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This  is  Volume  2  of  a  two  volume  tinal  repon  for  the  contract  "Analysis  of 
Dynamical  Plasma  Interactions  with  High  Vr)ltage  Spacecraft  "  The  period  of  technical 
performance  was  31  December  198X  through  31  Decetnber  1991,  The  objectives  of 
this  contract  were  to  study  dynamical  plasma  interactions  with  high  voltage  spacecraft, 
to  construct  a  three-dimensional  computer  code,  DynaPAC,  as  a  workbench  for  such 
studies,  and  to  support  the  SPEAR  program.  Volume  I  is  a  compilation  of  work  done 
to  model  high  voltage  plasma  interactions,  with  application  to  the  SPE.4K-1I  chamber 
tests  and  to  the  design  of  SPEAR  3.  Volume  2  contains  documentation  for  the 
DynaPAC  code,  as  well  as  for  the  two-dimensional  Gilbcn  code. 


vil 


1.  Introduction  and  Overview 


1.1.  Need  for  Dynal’AC 

At  the  inception  of  this  contract,  recent  flight  experitncnis  (notably  SPEAR-h  and 
concurrent  computer  modeling  demonstrated  a  strong  capability  to  predict  steady-state 
interactions  between  high  voltage  spacecraft  and  the  space  plasma.  However,  the 
time-dependent  response  of  the  space  plasma  to  the  high  lields  and  voltages  associated 
with  pulsed  power  systems,  and  the  associated  dynamic  spacecraft  charging,  had  not 
been  investigated.  Processes  not  adequately  modeled  included  formation  of  the  space 
charge  sheath,  current  flow  in  the  quasi-neutral  preshcath.  breakdown  phenomena, 
plasma  kinetics,  ionization  pi.K'esses.  and  the  effect  of  dynamic  processes  on  space 
craft  charging. 

This  inadequacy  became  apparent  in  trying  to  make  plasma  interaction  pre<licti(Mis 
for  the  SPEAR-ll  high  voltage  test  Equilibrium  sheath  calculations  gave  verv 
diflerent  results  from  sheath  calculations  using  shoil  time  approximations,  and  plasma 
currents  to  the  high  voltage  components  could  not  be  calculated  with  conlidence 
Thus,  it  was  decided  to  employ  S-CUBED's  experience  in  developing  equilibrium 
computer  models  to  develop  a  realistic  dynamic  computational  capability  ineoiporaling 
modern  and  advanced  computational  techniques. 

S-CUBEI)’s  three-dimensional  equilibrium  codes,  NA.SCAP/GEO. 
NAvSCAP/LEO.  and  POLAR,  have  proved  their  general  utility  through  application  to 
the  design  and  analysis  of  a  succession  of  laboratory  experiments,  spaceflight  experi¬ 
ments,  and  functional  spacecraft.  The  goal  of  DynaPAC  is  to  provide  a  similarly  gen¬ 
eral  tool  which  will  be  used  for  design  and  analysis  of  future  pulsed  power  systems  in 
space  and  dynamic  plasma  interactions  of  other  future  space  systems 

1.2.  DynaPAC  Capabilities 

DynaPAC  is  currently  a  user-friendly  and  programmer-friendly  workbench  for 
studying  plasma  interactions  with  realistic  spacecraft  in  three  dimensions.  Presently, 
DynaPAC  enables  plasma  interactions  specialists  to  perform  realistic  analyses  with 
direct  application  to  engineering  problems.  Its  current  capability  is  illustrated  by  its 
application  to  SPEAR-ll.  When  mature.  DynaPAC  will  he  usable  by  spacecraft 
engineers  with  plasma  interactions  expertise. 

The  core  capabilities  of  DynaPAC  are  to 

<  I )  Define  the  spacecraft  geometry  and  the  structure  of  the  computational  space: 

(2)  Solve  the  electrostatic  potential  about  the  (vbjcct.  with  flexible  boundary  condi 

lions  on  the  object  and  with  space  charge  computed  either  fully  by  partieles.  fully 

analytically,  or  in  a  hybrid  manner,  and 
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(3)  Generate,  track,  and  otherwise  process  representative  macroparticles  of  various 
species  in  the  computational  space. 

The  core  modules  are  designed  to  have  the  maximum  practical  user  control  and  to 
facilitate  the  incorporation  of  new  or  modified  algorithms.  Preprocessors  are  provided 
to  set  boundary  conditions  and  generate  input  tiles  in  a  modern,  screen-oriented  way. 
Similarly,  screen -oriented  postproce.ssors  are  provided  for  graphical  and  textual  data 
display. 


DynaPAC 

Software  Module  Structure 


Figure  1,1.  Modular  structure  of  the  DynaPAC  code. 
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1.3.  DynaPAC  Software  Klements 

Figure  1.1  shows  the  modular  structure  uf  the  DynaPAC  code  Except  for  the 
"Neutral  Gas  Model,"  all  modules  shown  are  currently  operatumal. 

DynaPAC  is  centered  around  a  DataBase  Manager,  The  DataBase  Manager  is  a 
library  of  routines  capable  of  making  large  arrays  of  information  contained  in  disk  fdes 
accessible  to  computational  modules.  It  has  a  programmer-friendly  language  for 
defining  data  types  and  for  retrieving  and  storing  data.  The  DataBase  Manager  will  be 
used  to  write  DynaPAC  results  in  the  format  required  for  other  application  codes,  and 
to  retrieve  information  generated  by  other  codes  for  u.se  in  DynaPAC  The  DataBase 
strategy  enables  DynaPAC  to  be  operable  on.  and  portable  among,  modem  high-power 
workstations,  which  have  proven  to  be  more  cost-effective  than  supercomputers  for 
this  type  of  code  development  and  analysis. 

.Spacecraft  geometrical  definitions  for  DynaPAC'  are  done  using  standard  finite 
element  pre-processors  such  as  PATRAN.  Among  the  advantages  of  this  approach  are 
that  the  geometry  can  be  realistically  represented,  and  that  finite  element  models  of  a 
spacecraft  constructed  for  other  purposes  can  be  adapted  for  DynaPAC  use.  The  com¬ 
putational  space  around  the  spacecraft  is  constructed  interactively  using  the  GridToo! 
module  Arbitrarily  nested  subdivision  allows  resolution  of  important  object  features 
while  including  a  large  amount  of  space  around  the  spacecraft.  The  DynaPAC  object 
definition  interface  program  PatDyn)  initializes  the  DataBase  and  con.structs  the 
properties  of  those  finite  elements  neighboring  the  spacecraft. 

DynaPAC  uses  a  high-order  finite  element  representation  for  the  electrostatic 
potential  that  assures  that  electric  fields  are  .strictly  continuous  throughout  space.  (The 
originally  proposed  tri-quadratic  finite  element  representation  proved  to  have  some 
unfortunate  properties,  and  was  replaced  by  the  new  scheme.)  The  electrostatic  poten¬ 
tial  solver  uses  a  conjugate  gradient  technique  to  solve  for  the  potentials  and  fields  on 
the  spacecraft  surface  and  through  the  sunounding  space  .Space  charge  options 
presently  built  in  include  Laplacian,  equilibrium  sheath,  "frozen  ions"  (appropriate  to 
the  early  aage  of  a  negative  transient  pulse),  "mobile  ions  -  barometric  electrons” 
(appropriate  to  the  several  microsecond  timescale  response  to  a  negative  pulse),  and 
"full  PIC"  A  screen-oriented,  menu-dnven  preprocessor  is  available  to  generate  initial 
conditions  and  to  build  an  input  file  for  the  potential  solver. 

Particle  tracking  is  used  to  study  sheath  currents,  to  study  particle  trajectories,  or 
to  generate  space  charge  evolution  for  dynamic  calculations.  DynaPAC  generates 
macroparticles  either  at  a  "sheath  boundary"  or  throughout  all  space.  Particles  are 
tracked  for  a  specified  amount  of  time,  with  the  timestep  automatically  subdivided  at 
each  step  of  each  particle  to  maintain  accuracy.  The  current  to  each  surlace  cell  of  (he 
spacecraft  is  recorded  for  further  processing.  A  future  task  is  to  tievelop  algorithms 
which  will  allow  DynaPA(  to  study  problems  having  neutral  ionization 
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A  screen -oiicntai,  tnenu-driven  pusiproccssui  is  used  U'  access  the  DataBase  Ides 
and  generate  graphical  output  illustrating  such  quantities  as  nhject  surlace  pt'ientiais, 
space  potentials,  particle  positions,  or  particle  trajcctcries.  Pdtential  contour  plots  may 
be  made  either  as  line  plots  or  as  color-filled  plois.  Contour  levels,  plot  range,  rn 
can  be  modified  through  the  user  interlace.  Plots  are  written  to  an  intermediate  graph¬ 
ics  file  which  can  be  displayed  using  any  of  several  terminal  protocols  t)r  converted  to 
PostScript  tiles  for  printing.  This  schente  is  designed  to  easdy  accommodate  new 
graphics  hardware  or  software. 


2.  DynaPAC  Users  Manual 


2.1.  Preparing  to  Use  DynaPAC 

This  manual  describes  the  use  of  DviiaPAC  Irorn  the  user  s  point  of  view. 
Because  DynaPAC  is  a  large  set  of  complex  codes,  tiiols  have  been  developed  to  make 
file  handling  and  input  generation  as  automated  an  well-interfaced  as  possible. 

DynaPAC  should  be  installed  on  a  reasonably  powerful  UNIX  workstation  The 
DynaPAC  program  tree  requires  about  50  MB  of  disk  space.  At  least  several  MB  of 
disk  space  should  be  available  for  problem  files.  Complex  problem  file  sets  can  easily 
re(juire  100  MB  or  more. 

2.1.1.  Units 

In  general.  DynaPAC  operates  in  the  SI  or  MKS  system  of  units  Electrostatic 
potentials  are  internally  stored  in  volts,  and  electric  fields  in  volts  per  meter.  Magnetic 
fields  are  always  in  tesla  (or  webers  per  square  meter)  Particle  energy  or  plasma  tem¬ 
perature  is  usually  in  electron-volts.  Charge  density  (p)  is  usually  divided  by  the  per¬ 
mittivity  of  free  space  so  that  it  has  the  units  of  volts  per  s(.|uare  meter. 

2.1.2.  DynaPAC  Environment 

DynaPAC  makes  use  of  UNIX  environment  variables  and  aliases.  These  are 
defined  by  the  "DynaPAC'. setup"  tile,  found  in  the  head  of  the  DynaPAC'  tree.  Figure 
2.1  shows  the  "DynaPAC. setup"  file  Only  the  first  three  "setenv"  commands  need  to 
be  customized  to  your  system;  the  .SCREENPKG  variable  needs  to  he  set  to  the  head 
directory  of  the  DynaPAC  screen  package,  the  DYNAPAC,'  vanable  needs  to  be  set  to 
the  head  di'cctory  of  DynaPAC'.  and  the  CPL'  TYPE  variable  needs  to  bo  set  to  either 
IRIX  (for  .Silicon  Graphics.  Inc.  computers  with  IRIX  4.0  or  later  operating  system)  or 
SLN4  (for  Sun  Microsystems  Sun-4  or  SparcStation  computers). 

DynaPAC  executables  are  suffixed  with  the  CPI!  TYPE  variable  Those  execut¬ 
ables  with  screen  interfaces  must  be  run  using  shell  scripts  which  set  additional 
environment  variables  prior  to  running  the  executable.  The  aliases  in  the 
DynaPAC'  setup"  file  make  all  this  transparent  to  the  user  It  is  not  necessary  (or 
recommended)  that  the  SDYN.APACYbin  directory  be  included  in  the  user's  "path".  It 
IS  recommended  that  the  DYNAPAC’  variable  be  defined  m  the  user's  '  .eshre"  file. 


n 

#  </;W'7r  'AC, 'A  s-crHi-;i) 

#  setup  new  dynapae  window 

#  checlt/update  the  shell  variables  below 

ff  shell  vanables 

#  DYNAPAC'  is  the  head  of  Dynapae  notle 

#  SCRPKNPKCi  is  the  head  of  ScreenPk*: 

#  CPfi  TYPl-l  is  iiiaehine  type 

U  '  (Sl'N4.  SI  NV  K«(K).  IRIS,  IRIX.  STI  I  I.IXl 
setenv  SCRhIiNPKti  /ileh/SereenPkg 
setenv  DYNAPAC  /iieh/dynapac 
setenv  CPI  TYPi  IRIX 

#  screen  package  things 

setenv  I  l-.RMC.AP  SS('KI',l  7NPK(i/src/teritk  ap/temicap 
setenv  f  ORMGIINDIR  VSCRIiPNPKCi/bin 


#  aliases 

alias  formgen  SFORMCiEiNDIR/formgen 
alias  Dmake  SDYNAPAC/dynamake 

#  executables 

aliits  Potent  SDYNAPA(7bin/Poieni_$CPL:  J  YPL 
alias  PatOyn  $DYNAPA(7bin/PatDyn_$rPl'  TYpp 
alias  PolDyn  $DYNAPAC/bin/PolDyn_$CPC jrYPli 
alias  NisDyn  $DYNAPAC/bin/NisDyn  SCPl  TYPE 
alias  PiiflGcn  $DYNAPA(7bin/Part(ien”sCPC  TYPE 
alias  Tnicker  $DYNAPA(7bin/Tracker  "SCPl  TYPE 
alias  ObjPotl  SDYNAPACVbin/ObjPotfsCPl  TYPE 
alias  MakeNpl  $DYNAPAC7binA1akcNpl_S(’Pl  _ I'YPE 
alias  AddGrid  SDYNAPA(7bin/AddGrid  Sf’PI  TYPE 

M  executables  with  screen  interface 
alias  DynaPost  $DYNAPA(7hin/DynaPost 
alias  DynaPre  $DYNAPAC/bin/DynaPre 
alias  GridTool  $DYNAPA(7bin/Gridl  ool 
alias  Scanner  $DYNAPA(7bin/Scanner 

#  drivers 

alias  dbtool  SDYNAPA(7src/dblib/dblool 

#  greetings 
echo  ' 

echo  Welcome  to  DynaPAC  . 
echo 

cd  SDYNAPAC 
pwd 

Figure  2.1.  "DynaPAC  setup"  tile 

Tt)  set  up  the  aliases  and  envinniineiit  variables  needed  to  run  DynaPAC.  change 
directories  to  the  DynaPAC  head  directory  and  type 
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source  DynaPAC. setup 


or  (from  any  directory) 

source  SDYNAPAC/DynaPAC. setup 
You  will  be  left  in  the  DynaPAC  head  directory 

2,1.3.  DynaPAC  Modules 

DynaPAC  is  a  collection  of  executable  modules  which  perform  the  core  functions 
of  DynaPAC  and  associated  pre-  and  post-processing.  The  first  module  that  a  user 
will  encounter  is  not  part  of  DynaPAC  at  all,  but  an  external  finite  element  geometry 
generator  (presently  either  PATRAN  or  POLAR)  used  to  define  a  geometrical  object 
for  DynaPAC  simulation.  Then  DynaPAC  modules  will  be  applied  to  the  object  in 
more  or  less  the  following  order; 

GridTool  is  used  to  interactively  define  an  arbitrarily  nested  cubic  grid  system 
for  the  space  surrounding  the  object. 

DynaPre  is  used  to  ( I )  define  input  for  the  geometry  interface  modules.  PatDyn 
and  PolDyn.  DynaPre  will  later  be  revisited  to  (2)  define  initial  potential  boun¬ 
dary  conditions  on  the  object  surfaces;  (3)  define  input  to  the  potential  solver 
module.  Potent;  (4)  define  input  to  the  particle  generator  module.  PartGen;  and 
(5)  define  input  to  the  particle  tracker  module.  Tracker. 

PatDyn  reads  a  PATRAN  neutral  file  and  generates  the  surface  and  volume  infor¬ 
mation  required  by  the  DynaPAC  computational  modules. 

PolDyn  reads  geometrical  information  from  POLAR  program  files  and  generates 
the  surface  and  volume  information  required  by  the  DynaPAC  computational 
modules. 

Potent  is  used  to  calculate  the  electrostatic  potential  field  in  the  space  surround¬ 
ing  the  object. 

Scanner  is  primarily  used  to  display  the  electrostatic  potential  field  in  the  space 
surrounding  the  object  It  can  also  be  used  to  display  selected  other  spatial  vari¬ 
ables.  and  to  print  any  valid  DynaPAC  spatial  or  surface  variable. 

Part(»cn  is  used  to  generate  macroparticles  to  study  particle  trajectories,  surface 
currents,  and  space  charge. 

Tracker  is  used  to  track  macropart  ides  to  study  particle  trajectories,  surface 
currents,  and  space  charge, 

DynaPost  is  a  postprocessor  to  display  (and  plot  time  history  of)  values  of  sur 
face  currents,  and  to  generate  input  for  (and  run)  the  ObjPotl  module, 

ObjPotl  plots  isometric  views  of  the  object  surface,  color-coded  by  material 
number,  conductor  number,  surface  potential,  surface  electric  field,  or  incident 
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flux. 


2.1.4.  DynaFAl'  Files 

The  user  brings  to  DynaPAC’s  ClritlTool  module  files  describing  the  geomeinc 
and  material  configuration  of  his  or  her  object.  At  that  time  a  prefix  will  be  assigned 
to  the  file  set.  From  that  point,  all  files  are  created  by  the  DynaPAC  modules.  In  gen¬ 
eral.  the  following  naming  conventions  arc  followed: 

prefix.suffix  files  (where  suffix  is  not  capitalized)  are  CridTool  ascii  output  files. 
They  may  be  edited  by  the  user  (though  altering  them  is  not  recommended) 

prcfix.SUFFIX  files  (where  SUFFIX  is  capitalized)  are  DataBase  Manager  files 
In  general,  these  are  non-aseii  and  cannot  be  user  iiiodified  E^ach  such  file  has  a 
maximum  size  of  MB, 

name  in  tiles  (where  name  identifies  a  module)  are  ascii  input  files  They  con¬ 
tain  well-labeled  parameters  whose  values  may  be  altered  by  the  user. 

name  out  files  (where  name  identifies  a  module)  are  output  files  These  consist 
primarily  of  reports  of  the  progress  of  the  module,  but  also  contain  information  of 
interest  to  the  analyst.  Various  modules  also  write  auxiliary  output  files 

Scratch.XX  files  are  created  by  the  Potent  module  and  may  be  deleted  any  time 
Potent  is  not  operating. 

fort.2  files  are  graphics  files,  and  are  overwritten  each  lime  a  figure  is  generated 
by  a  DynaPAC  module.  Figures  can  be  saved  by  renaming  the  corresponding 
fort.2  file. 

More  specifically,  here  is  a  list  of  some  important  DynaPAC  tiles  with  their  con¬ 
tents: 

prefix.neu  is  a  PATRAN  neutral  file  (ascii)  containing  object  geometrical  infor¬ 
mation. 

prefix.mat  is  a  user-generated  file  (ascii)  containing  material  definitions. 

prefix.grd  (ascii)  is  generated  by  CridTool  to  define  the  DynaPAC  grid  system 
for  the  problem. 

prefix.obj  (ascii)  is  generated  by  CridTool  as  a  simplified  object  description. 
prefix.DP  is  the  main  DataBase  file  for  the  problem. 

prefix.HI  (ascii)  contains  the  data  packet  names  and  definitions  for  the  DataBase 
Manager. 

prefix.BS  and  prcfix.MKnn  are  generated  by  PatDyn  or  PolDyn  to  contain  spe¬ 
cial"  element  information.  These  files  are  treated  as  read-only  by  subsequent 
modules,  and  can  be  shared  among  tlifferent  versions  of  a  problem  using  the  same 
geometry  and  gridding.  Also,  the  prefix.MEn[i  files  are  read  only  by  the  Potent 
module,  and  may  be-  deleted  if  no  further  potential  solutions  are  anticipated. 
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prefix.PTn  files  arc  used  to  slorc  currcnj  nuica)paniclc  inlDrinaiion. 

prefix. POTG.  prefix-I^OTS,  and  prefix.GlJR  store  time  histories  of  spatial  poten¬ 
tials,  surface  potentials,  and  surface  currents. 

User-modifiable  input  files  include 


Default  Name 

Module 

Generated  By  I 

caddy n  in 

PatDyn 

DynaPre  j 

caddyn  in 

PolDyn 

DynaPre  | 

ps  in 

Potent 

DynaPre  I 

pgjn 

PartGen 

DynaPre  j 

tr  in 

Tracker 

DynaPre  ! 

tr  traj  in 

Tracker 

DynaPre 

1  obipull  in 

ObjPoll 

DynaPosi 

.....  i 

.Names  for  standard  output  tiles  are  defined  by  the  user. 
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2.2.  Handling  Input  and  Output 

The  modules  of  DynaPAC  are  conveniently  divided  into  computational  module', 
(such  as  PatDyn.  Potent,  PartGen,  Tracker,  and  ObjPotl).  and  interactive  modules  (such 
as  GridTool.  DynaPre,  Scanner,  and  DynaPost),  The  computational  modules  read  ihcir 
parameters  and  directives  from  standard  input  and  write  (mostly  diagnostic)  informa¬ 
tion  to  standard  output.  For  these  modules  it  is  convenient  to  obtain  an  initial  input 
stream  (from  DynaPre  or  DynaPost)  and  use  the  I’NIX  editor  for  subsequent 
modificatio.is  ol  parameter  values.  The  interactive  modules  obtain  user  input  through 
the  screen  handler,  and  produce  output  via  screen  display,  graphics  hies,  or  ascii  tiles 
We  assume  the  user  is  familiar  with  the  standard  UNIX  utilities,  and  describe  here 
how  to  use  the  DynaPAC  Screen  Handler  and  the  graphics  display  interfaces. 

2.2.1.  Using  the  Screen  Handler 

The  screen  handler  provides  a  method  of  displaying  parameter  values  to  the  user 
in  a  convenient  and  well  organized  form,  while  allowing  him  or  her  to  alter  the  \alues 
directly  on  the  form.  The  scheme  wurks  on  terminals,  terminal  windows,  or  terminal 
emulators 

2.2.1.1.  Appearance  of  the  Screen 

Figure  2.2  shows  a  typical  screen  On  the  top  line  appears  the  static  top  level 
menu  for  the  module  being  run.  Below  that  (in  reverse  video  where  available)  is  the 
"annunciator  line",  which  occasionally  contains  useful  messages  about  the  process 
being  performed.  Second  and  third  level  menus  and  screens  are  superimposed  on  one 
another,  with  only  the  lowest  level  menu  accessible  to  the  user.  Menus  and  screens 
vary  in  size  from  small  pulldowns  to  the  entire  area  below  the  annunciator  bar. 
depending  on  the  space  required  for  their  function. 


winterm 
Exit  Print 
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Print 

WORKING  .... 

PLOT 

EDITOR 

Output  Mode: 

FILE 

Data 

File  Prefix: 

slabs 

Output  File  Name: 

Scannet  . out 

Data 

File  Type: 

DYNAPAC 

Data  Representation:  32_NODE 

Data 

None; 

POT_Grid 

Hiixlou  Manager: 

xn 

Data 

Type: 

SPATIAL 

Plotter : 

GLDr on 

Cut  Plane  offset:  5.0000  Prinary  Grid:  Mx=  *)  My=  9  N2=  9 

Diaq  Level;  1  Cut  Plane  Normal:  Y 

Horizontal  Axis;  X 

-  Options  -  Vertical  Axis:  Z 

GRID  LIMITS  LABELS  i  INCLUDES 


Figure  2.2.  Example  of  DynaP AC  screen. 


2.2. 1.2.  Types  of  Fields  (with  Examples) 

Several  types  of  tields  appear  on  the  screens 

Static  Fields  are  used  to  describe  or  iiluininatc  the  ct)ntenis  of  the  screen,  and  are 
not  selectable  by  the  user.  For  example,  the  ntaximum  grid  dimensions  arc 
displayed  in  the  "Plot  Edit  Grid  Limits"  submenu  for  user  convenience, 
and  may  not  be  changed. 

Value  Fields  contain  numerical  or  ascii  values.  The  values  may  be  modified  by 
positioning  the  highlight  over  the  value  field,  selecting  the  field  (with  a 
<CR>),  typing  the  new  value,  and  affirming  the  entry  (with  a  second  <CR>) 
The  screen  handler  performs  some  validity  checking  on  value  field  entries 
Examples  of  value  fields  are  "Plasma  Density''  and  "Plot  Title  ' 

Toggle  Fields  are  value  tields  with  a  small  number  of  predefined  valid  entries 
Toggle  fields  are  modified  by  selecting  the  field  (as  above),  using  the  spa 
cebar  to  "toggle"  through  the  possible  options,  and  selecting  the  desired 
option  (with  a  <CR>)  An  example  of  a  toggle  field  is  "UNIT.S",  which  may 
lake  the  values  "METER.S".  "CENTIMETERS".  "IN'CHES  .  FEET  ,  or 
"OTHER".  (In  this  case,  choosing  "OTHER"  results  in  display  of  a  value 
field  to  enter  (he  desired  unit  in  meters.) 

Menu  Fields  result  in  display  of  a  lower  level  menu.  Selecting  the  "Edit  Scnpi ' 
field  results  in  display  of  the  input  parameters  for  a  computational  module 
for  user  inspection  and  modification.  Similarly,  selecting  the  ’Edit  Plot" 
field  allows  the  user  to  choose  the  parameters  of  a  plot. 

Action  Fields  result  in  some  action  taking  place.  "Read"  causes  previous 
relevant  parameter  or  option  values  to  be  read  from  DynaPAC  files,  'Write" 
causes  the  results  of  an  interactive  session  to  be  written.  "Make  Script'  or 
"Make  Plot"  cause  an  input  stream  or  plot  to  be  generated  and  written  on  the 
appropriate  file  "Show  Plot "  causes  spawning  of  a  shell  to  run  the  currently 
selected  graphical  interface  program.  "Exit'  results  in  escape  to  the  parent 
menu,  and  "Help"  sometimes  gets  display  of  a  rudimentary  help  file.  Aefion 
fields  of  the  form  "Run  Program"  generally  do  not  work  and  should  be 
avoided. 

2.2. 1. 3.  Navigating  the  Screen  and  Menu  Tree 

The  screens  und  the  menu  tree  arc  navigated  using  the  letter  keys  "h".  'j'  ,  "k  ’. 
and  "I",  the  Einter  <(’K>  key.  and  the  l-Acape  <E:SC>  key 

The  letter  keys  move  (he  highlight  respectively  left,  down,  up.  and  right  on  the 
current  screen  or  menu.  In  general,  if  (he  requested  molii'n  is  not  possible  the  screen 
handler  will  perform  an  alternate  motion  (e  g.,  up  instead  of  left)  or  will  wrap  around. 


It 


The  Enter  <CR>  key  is  used  to  select  a  field  or  to  affirm  the  niodified  value  of  a 
previously  selected  field. 

The  Escape  <ESC>  key  is  used  to  rest<»re  the  previous  value  of  a  currently 
.selected  field,  or  (if  no  field  is  selected)  to  signify  con)pletion  of  action  on  the  current 
menu  or  screen.  Escape  from  a  menu  is  usually  to  the  parent  menu,  but  may  result  in 
display  of  another  menu  at  the  same  or  a  lower  level  if  the  context  is  appropriate 
Some  menus  will  not  honor  the  <ESC>  key.  but  require  selection  of  the  "Exit"  field. 
In  particular,  the  top  level  menu  does  not  accept  the  <ESC>  request. 


Key 

h 

.) 

k 

1 

<CR> 
<CR> 
<ESC> 
I  <ESC> 


Action 

Move  left  (or  up) 

Move  down  (or  right) 

Move  up  (or  left) 

Move  right  (or  down) 

Select  field  (for  modification  or  action) 
Accept  new  value  for  selected  field 
Restore  old  value  for  selected  f.  .d 
Signify  menu  completion 


2.2.2.  Displaying  Graphics 

DynaPAC  modules  write  graphics  commands  to  the  file  "fort. 2".  from  which  they 
may  be  read  and  executed  by  any  of  several  graphical  interface  programs.  The  graphi¬ 
cal  interface  programs  may  be  executed  by  the  "Show  Plot"  menu  selection  from  the 
interactive  modules,  or  directly  from  the  DNIX  shell.  In  the  former  case,  the 
correspondence  between  plot  interface  name  and  executable  interface  is  defined  by  the 
SDYNAPAC/bin/define_plot_readers  file.  Note  that  the  "Make  Plot"  selection  appends 
a  plot  to  the  "fort. 2"  file,  and  the  "Show  Plot"  selection  resets  the  file  pointer  to  the 
beginning  of  the  file.  Thus,  the  first  "Make  Plot"  selection  following  a  "Show  Plot" 
selection  elTccliveiy  erases  all  previous  plots 

Currently  favored  graphical  interface  programs  are  listed  below.  Where  not  oth¬ 
erwise  specified,  frame  advance  is  accomplished  by  a  <CR>. 

GLDraw  or  IRIX  I’lot  displays  plots  through  calls  to  the  SGI  Gl,  library.  The 
plot  window  may  be  moved  or  resized  (with  rescaling  taking  effect  on  the 
following  frame)  Erame  advance  uses  the  right  mouse  button 

XDraw  displays  plots  through  calls  to  the  XLib  and  XToolkit.  The  pkit  window 
may  be  moved  or  resized  (with  rescaling  taking  effect  on  the  following 
frame).  Frame  advance  uses  any  mouse  button 

Tek4xxx  (where  4\\x  is  4014,  4105.  4107.  4207)  sends  to  standard  output  graph 
ics  eomiiiaiuls  appropriate  to  the  s|H‘eitied  Tektronix  terminal  The  frame 
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advance  coriimand  is  <CR>. 

PostGray  and  PostColor  generate  PostScript  commands  on  standard  output  Hie 
former  replaces  the  color  scale  for  contour  levels  with  a  monotonie  gra\ 
scale.  (No  pause  for  frame  advance.) 

Adobe  generates  (on  stamlard  output)  commands  which  can  be  iviade  acceptable 
for  input  to  the  Adobe  Illustrator  program  (No  pause  for  frame  advance  ) 

SunCore  and  SunPix  display  graphics  commands  in  the  SunView  windowing 
environment.  The  frame  advance  command  is  <rR>. 


i.S 


2.3.  Defining  Objects  and  (>rids 

The  first  two  tasks  of  a  DynaPAC  simulation  are  to  define  a  spacecraft  or  test 
object  as  a  boundary  surface  representation,  and  to  define  a  grid  structure  about  the 
object.  DynaPAC  provides  no  object  definition  facility,  but  accepts  definitions  pro¬ 
duced  by  other  codes.  Presently  supported  are  PATRAN  (a  commercially  available 
finite  element  analysis  code)  and  POLAR  (a  spacecraft  charging  analysis  code 
developed  by  S-CUBED  fur  the  Air  Force).  Other  geometry  definition  codes  may  be 
supported  as  requirements  dictate. 

DynaPAC’s  GridTool  module  is  used  to  define  an  arbitrarily  nested  grid  structure 
about  the  object.  With  the  definitions  of  object  and  grid  established,  the  PatDyn  (for 
PATRAN  objects)  or  PolDyn  (for  POLAR  objects)  module  is  run  to  analyze  the 
geometrical  configuration  and  establish  the  DynaPAC  database. 

2.3.1.  Defining  Objects  u.sing  PATRAN 

Object  definition  for  DynaPAC  is  identical  to  object  definition  for  NASCAP/LEO. 
For  tutorial  exercises  on  using  PATRAN  to  define  objects,  the  reader  is  referred  to  the 
User's  Guide  to  NASCAPll.lX).  Here  we  give  only  the  specific  requirements  that  the 
object  must  satisfy. 

The  user  must  construct  the  object  surface  using  linear  triangles  (TRI/.^/icond) 
and  linear  quadrilaterals  (QU AD/4/icond).  The  PATRAN  "configuration  code"  is 
interpreted  as  the  "conductor  number";  the  PATRAN  MID  is  used  as  the  material 
number.  The  object  should  have  no  free  boundaries,  and  all  surface  normals  should  be 
consistently  outward.  Nodes  should  be  equivalenced  as  appropriate,  and  "COMPACT" 
and  "ELEM  COMPACT"  operations  .should  he  performed  prior  to  writing  the  neutral 
tile. 

The  PATRAN  neutral  file  (output  as  "patran.out./t")  should  be  renamed  to 
prejix.neu" .  A  material  file  which  a.ssociates  material  names  and  properties  with  the 
material  numbers  assigned  using  PATRAN  must  appear  as  "prefix.mal" .  The  material 
file  contains  two  types  of  packets: 

(1)  A  one-line  packet  consisting  of  an  integer  (MID)  followed  by  a  material  name. 
PatDyn  then  associates  the  material  name  with  all  surfaces  have  the  MID.  If  the 
material  name  is  defined  in  PatDyn’s  material  database,  its  default  properties  are 
assigned.  If  not.  a  warning  is  printed. 

(2)  A  four-line  packet  consisting  of  a  material  name  (which  must  have  already  been 
assigned  an  MID)  followed  by  three  lines  of  materia!  properties  in  NASCAP/LEO 
format. 

Figure  2.3  shows  an  example  material  file. 

1  Alum 

2  Kapton 
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3  Teflon 

4  NPaint 

5  Silver 

6  CPaint 

7  Gold 

8  Aquadg 
Gold 

1  .(X),.00 1 1  ..79...88..8,88.79.,92 

53.48, 1 .73,4 1 3, 1 35..-00(K)29.- 1 ..  I  .E+4.2.E+3 

l.E-13,l..l.E+3,2(). 

Figure  2.3.  Example  material  file,  using  8  pre-delined  material  names,  then 
redefining  the  properties  of  "Gold." 

2.3.2.  Defining  Objects  using  POLAR 

It  is  simpler,  although  more  restrictive,  to  use  POLAR  to  define  DynaPAC 
objects.  The  procedure  is: 

(1)  Define  a  valid  POLAR  object  using  POLAR’s  VEHICL  module.  Be  sure  that  ail 
POLAR  checks  are  passed.  Be  sure  to  correctly  define  the  DXMESH  parameter 

(2)  Keep  POLAR  output  files  "fort. 21",  "fort.  17",  and  "fort.  1 1”.  Rename  "fort. 21”  to 
"fort.8"  (mv  fort.2l  fort. 8).  Make  extra  links  to  these  files  to  avoid  their  being 
deleted. 

(3)  No  extra  material  file  is  needed.  The  material  names  and  definitions  will  be  read 
directly  from  the  POLAR  files. 

(4)  Run  GridTooI  and  PolDyn  as  described  below. 

2JJ.  Creating  Grids  using  (iridTtMtl 
2J3.1.  General 

GridTooI  allows  the  user  to  define  a  grid  structure  for  an  object  created  by 
PATRAN  or  POLAR.  It  runs  in  a  "screen  oriented  input"  mode  on  UNIX  machines 
with  standard  A.Sf’ll  terminals  or  emulators,  and  displays  graphics  using  any  of  the 
plot  di.splay  interface  modules.  A  .Silicon  Graphics  Iris  version  with  interactive  graph¬ 
ics  was  developed  for  an  old  SGI  configuration,  but  has  not  been  updated  to  the 
current  SGI  windowing  system. 

Many  of  the  forms  or  screens  shown  in  this  file  are  slightly  out  of  date.  and.  in 
particular,  the  graphical  options  do  not  app'^ar.  The  differences  in  the  current  forms 
should  be  self  explanatory. 


17 


2  JJ.2.  I/O  Files 

GridToo!  reads  the  neutral  file  (suffix  ncu)  (renamed  from  "patran.out.n")  from 
PATRAN,  or  the  files  fort.l  I.  fort.  17,  and  foil  S  (renamed  from  fo!l,2l)  from  POLAR 
It  creates  two  output  tiles  (suffices  .grd  and  .obj).  The  output  files  are  also  read  in  for 
a  restart  run. 


File 

Contents 

<pref!X>.ncu 

neutral  file  from  PATRAN 

fort.  1 1 

POLAR  file 

fort.  1 7 

POLAR  file 

fort.  8 

POLAR  file  fort.2l 

<pretix>.grd 

grid  parameters  (ASCII  text) 

<prefix>.obj 

surface  information  (ASCII  text) 

The  user  will  specify  <prefix>  which  will  be  the  same  for  all  three  files.  The 
prefix. grd  file  will  be  read  by  PatDyn  or  PolDyn  to  define  the  grid  structure.  The 
prefix. obj  file  is  used  only  for  restarts  of  GridTool. 

2.3.3.3.  GridTool  Capabilities 

GridTool  presently  has  the  following  capabilities; 

1.  Create  and  modify  a  primary  grid  around  a  PATRAN  or  POLAR  object. 

2.  Add  and  modify  a  child  grid. 

Delete  a  grid  along  with  all  of  its  child  grids. 

4.  Restart  from  an  existing  grid  structure. 

5.  Graphically  or  tubularly  display  the  current  grid  structure. 

2.3.3.4.  Running  GridTool 

It  is  best  to  illustrate  how  to  use  GridTool  by  going  through  an  example.  Let's 
say  we  have  a  neutral  file  generated  by  PATRAN,  and  the  file  name  is  Cube.neu' 
Now.  start  up  GridTool.  The  first  screen  will  be  displayed  with  the  following  menu: 


I  Exit 

Files 

Edit 

Show 

Write 

_ 

Reset 

Hei^ 

IX 


Should  the  Help  field  be  selected,  the  tbllowing  screen  will  be  displayed; 


j  Grid  fool  Help  ij 

; 

1 

Exit 

Leaving  GridTool  p 

Files 

Set  file  prefix 

Show 

Show  the  grids  dclined  so  far  i 

Edit 

Edit  grids  structure  | 

Show 

show  parameters  of  current  grid  | 

Add 

add  new  child  grid  to  the  current  grid 

Delclc 

delete  current  grid  and  its  children  (if  any)  ; 

Modify 

modify  paiametcrs  of  current  grid  i 

Write 

Write  current  grids  to  file  <Prefix>.grd  j 

Reset 

Reset  everything  to  defaults  (to  start  all  over)  ; 

[H-lp 

fhis  screen  | 

Normally,  the  sequence  of  commands  selected  are:  (I)  Files;  (2)  Edit;  (3)  Show. 
(4)  Write;  (.S)  Exit.  Let's  now  go  through  these  steps. 

2.3J.4.1.  Flle.s 

Select  Files'.  The  default  prefix,  fort'  is  shown.  Hit  <KET>  to  change  it,  and 
type  in  the  correct  file  prefix  which,  in  this  case,  is  'Cube  .  1  hen  hit  <ESC>  to  exit 
this  menu.  A  new  form  with  two  toggle  ficKIs  is  then  popped  up: 


Neutral  file  type: 

FATRAN 

Units  used: 

meters  j 

Let's  say  the  units  used  to  define  the  object  was  in  inches.  We  then  would  move  to 
field  saying  meters'  and  to  change  it  to  'inches'.  Again,  hit  <ESC>  to  exit  this  torn. 

2.3J.4.2.  Edit 

Select  Edit'.  The  default  current  grid  (parent  grid)  is  grid  #1.  Since  wc  have 
not  defined  any  grid  yet,  start  from  grid  #1  is  a  must  Hit  <ESC>  to  exit  this  form  Wc 
will  then  be  in  the  following  pop-up  menu: 
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Exit  i 
Show  j 
Add  i 
Delete  ' 

Modify  I 

Here 

Exit  will  exit  this  menu  (^ame  as  hitting  <ESC>)  j 

Show  will  show  eurrent  grid  parameters  j 

Add  will  add  a  new  child  grid  to  current  grid  ! 

Delete  will  delete  the  <'urrent  grid  and  its  children  | 

(cannot  delete  primary  grid)  | 

Modify  will  modify  the  current  grid  | 

(can  only  modify  empty  (childless)  grids)  | 

Let's  select  Show  ui  see  the  current  grid  parameters.  The  following  form  will 
be  shown: 


I 

Current  grid  # 

Mesh  Size  (m) 

0.0254 

! 

Dimension 

5  5  5 

1  Child  Gnd  U  | 

Hit  <ESC>  to  exit  this  form,  then  select  'Modify  to  make  some  changes.  We  will  see 
the  following  Modify  form: 


2U 


We  now  can  move  around  the  form  making  changes.  Say  we  want  to  change 
grid  dimensions  to  10x10x20,  mesh  size  to  20cm.  The  form  now  shows  new  extents 
of  the  object  within  the  grid.  Let's  also  say  we  want  to  translate  the  object  8  mesh 
units  in  the  -Z  direction.  Select  TRANSLATE  OBJECT',  change  Dz  to  -8,  hit 
<ESC>.  Finally,  Select  'ACCEPT'  CHANGES'  to  save  the  changes  made  (or  select 
'IGNORE  CHANGES'  to  restore  previous  parameters  of  this  grid.) 

Enough  for  modifying  grids,  now  let's  add  a  child  grid.  Select  'ADD'  to  see  the 
ADD  form: 
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New  child  grid  # 

2 

Parent  grid  dimension; 

10  10  20 

Grid  limits  on; 

(Parent  Grid) 

(Primary  Grid) 

5  <=  X  <=  7 

5.0000  <=X<=  7.0(K)0 

5  <=  Y  <=  7 

5.0000  <=  Y  <=  7.00{X) 

10<=Z<=  12 

lO.fXXX)  <=  Z  <=  12.0(XX) 

Subdivision  Ratio: 

2 

RECALCULATE 

CHECK  OVERLAPS 

ACCEPT  GRID 

IGNORE  GRID 

. J 

A  set  of  default  parameters  have  been  set  to  this  newly  created  child  grid.  Some 
parameters  of  parent  grid  are  induced  in  the  form  to  help  us  figuring  out  the  location 
of  the  child  grid.  Again,  we  can  move  around  the  form  making  changes.  Other  func¬ 
tions  include; 


RECALCULATE 

update  the  form. 

CHECK  OVERLAPS 

check  overlaps  betw'een  child  grids 

(not  relevant  here  because  we  only  have  1  child  grid). 

ACCEPT  GRID 

save  this  child  grid 

IGNORE  GRID 

ignore  this  child  grid. 

Let's  say  we  are  satisfied  with  2  grids.  Select  'Kxii'  to  exit. 

2.3.3.4.3,  Show 

The  SHOW  option  allows  either  tabular  or  graphical  display  of  the  currently 
defined  grids.  In  this  case,  we  will  illustrate  the  TABLE  option. 

Select  'Show'.  Then  select  TABLE  from  the  pulldown  menu.  The  default  is  to 
show  parameters  of  all  defined  grids,  but  we  may  choose  to  look  at  some  subset  as 
well.  Exit  current  form  to  see  the  defined  grids: 
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There  are  2  grids  defined  j 

Grid  # 

1 

2 

Nx 

10 

7 

Ny 

10 

7 

Nz 

20 

7 

Parent 

0 

1 

Mesh  Ratio 

1 

3 

Mesh  Size 

0.02000 

Limits(in  parent  grid) 

Range  X 

1.  10 

5,  7 

Range  Y 

1,  10 

5,  7 

Range  Z 

1.  20 

mmm\ 

We  could  go  back  to  'Edit'  again  to  make  changes  if  we  are  not  pleased  with  the  grids 
shown.  Otherwise,  select  'Write'  to  write  new  grid  structure  to  file.  This  will  result 
in  2  files  being  created:  Cube.grd  and  Cube.obj. 


2.3.3.4.4.  Exit 

Select  'Exit'  to  exit  GridTool.  Now,  we  have  three  files: 

Cube.neu  neutral  file  (input) 

Cube.grd  grid  parameters  file  (restart  file) 

Cube.obj  surfaces  data  file  (resian  file) 

The  Cube.grd  will  be: 

2 

)  10  10  20  0  1  10  1  10  1  20 

1.00000  10.00000  1.00000  10.00000  I.IKXKX)  20.(X)no  2.00000E02 

1  1 

0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0  0  0.  0. 

2  7  7  7  1  5  7  .S  7  10  12 

5.000(K)  7.00000  .S.OIKKXI  7,tX)fl00  lO.tKKXX)  12(XMX1  6.66667E4)3 

3  3 

0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0.  0  0.  0  0.  0.  0.  0.  0. 

2.54000E  02  2,54(XX)E-02  2.-S4000E-02  0.  0.  O.IWXXM)  2..34000E-02 


23 


The  interpretation  of  the  prefix. grd  tile  is; 


Linefs) 

Fteld(s) 

Contents 

1 

Total  number  of  grids  defined 

2-5 

Grid  1  parameters 

2 

1 

Grid  Number 

2 

2-4 

nx.  ny.  nz 

2 

5 

parent  grid  number 

2 

6-!l 

limits  in  PARENT  grid  coordinates 

3 

1-6 

limits  in  PRIMARY  grid  coordinates 

3 

7 

mesh  size  ( meters  | 

4 

1 

mesh  ratio  with  respect  to  PARENT  grid 

4 

2 

mesh  ratio  with  respect  to  PRIMARY  grid 

5 

reserved  for  future  use 

6-9 

Grid  2  parameters  (as  2-5) 

10 

Object  size  in  meters 

10 

Object  center  relative  to  primary  grid  center  [meters] 

10 

mm 

Unit  conversion  factor 

2.3.3.5.  Afterthought  Object  Translation 

A  common  occurrence  is  that,  after  a  grid  structure  is  mostly  or  fully  defined,  the 
user  will  decide  that  a  grid  translation  is  needed.  This  can  be  done  in  the  files 
prefix. grd  and  prefix. obj: 

1.  Write  out  the  grid  structure  and  exit  GndTool. 

2.  Modify  fields  4-6  of  the  last  line  of  the  prefix. grd  tile  to  refiect  the  desired  object 
translation  fmeters|. 

3.  Modify  fields  4-6  of  the  second  line  of  the  prefix. obj  file  exactly  as  in  step  (2). 

4.  Restart  GridTool  to  display  and  check  the  new  object  position. 

2.3.4,  Processing  Objects  with  PatDyn  or  PolDyn 

PatDyn  or  PolDyn  should  be  run  in  a  directory  containing  only  the  output  files 
from  the  object  definition  program  and  from  GridTool.  plus  a  small  input  file  (which 
can  be  generated  with  the  CadDyn  option  of  the  DynaPre  module).  Normally,  only  the 
PREFIX  card  is  needed.  I'he  other  options  are  for  diagnostic  purposes. 

The  input  options  to  PatDyn  and  PolDyn  are: 

COMMENT 

Ignore  this  card. 

DIAGNOSTICS 

Print  (a  large  amount  of)  diagnostic  output. 
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ECHO 

EcIk)  the  neutral  file 

END 

Conclude  option  input. 

HELP 

Print  the  availahle  options 

NOECHO 

Do  not  echo  neutral  file  (I  his  is  the  default  ) 

NOPR(X'[-S.S 

Stop  before  generating  any  matrix  elements 

OFFDIAC. 

Suppress  off-diagonal  elements  between  surfaces 

PREFIX  prefix 

Deline  the  tile  prefix  for  this  run. 

PROCESS  IX  iy  k/.  GRID  ignd 

Proceed  as  expeditiously  as  possible  to  process  the  requested  special  clement 

(for  diagnostic  purposes). 

REMARK 

Ignore  this  card. 

PatDyn  or  PolDyn  will  write  an  output  file  of  diagnostic  information  which  you 
may  use  to  monitor  its  progress.  It  will  initiali/.e  the  DataBase  files,  and  create  the 
prefix. DP,  prefix. HI,  prefix. BS,  and  prefix. MEnn  (one  for  each  grid  with  special  ele¬ 
ments)  files. 

Early  in  execution  of  the  object  definition  interface  tile,  an  important  quantity  to 
monitor  is  "N.Spec  "  (You  may  "grep"  for  it.)  This  is  the  initially  determined  number 
of  special  volume  elements.  The  current  estimate  is  printed  three  times  for  each  gnd. 
The  limit  on  special  volume  elements  is  4(W5 

Many  of  the  "error  encountered"  messages  correspond  to  cells  originally  thought 
to  be  special,  but  on  careful  examination  turned  out  to  be  full  or  empty.  This  message 
is  particularly  likely  when  an  object  surface  is  coincident  with  a  grid  plane  The  final 
determination  (of  special,  tilled,  or  empty)  is  usually  correct. 


2.4.  Calculating  Potentials 

The  Potent  module  of  DynaPAC  is  used  to  calculate  potentials  in  the  space 
around  the  spacecraft.  Normally,  all  surface  potentials  are  held  at  fixed  potential 
However,  it  is  possible  to  specify  fixed  electric  field  at  insulating  surfaces,  or  to  allow 
spacecraft  potentials  to  change  in  response  to  plasma  currents.  .Several  ways  of  deal¬ 
ing  with  space  charge  are  provided,  and  advanced  users  are  encouraged  to  extend  these 
capabilities. 

2.4.1.  Using  DynaPre  (IPSl  to  Specify  Initial  Potentials 

When  the  IPS  option  is  chosen  from  the  DynaPre  top  menu  bar,  the  following 
options  are  presented: 

Fixit 

Return  to  the  DynaPre  top  menu  bar 

Read 

Read  initial  potentials  previously  written  to  the  DynaPAC  DataBase  files. 
Write 

Write  initial  potentials  into  the  DynaPAC  DataBase  files  This  option  must 
be  chosen  prior  to  exiling  IPS.  or  the  new  specification  will  not  take  effect. 

Conductors 

Specify  initial  potentials  for  conductors 
insulators 

Specify  initial  potentials  for  insulators 
Show  Table 

Show  the  initial  potential  table. 

Reset 

.Set  all  cells  to  fixed  potentials  of  zero  volts. 

Help 

No  help  is  currently  implemented  b>r  IPS 

The  most  common  initial  condition  setting  is  to  set  all  conductors  to  fixed  poten¬ 
tials  This  is  done  with  IPS  as  follows: 

(  I  I  Choose  Conductors”  from  the  pull-down  menu  Select  with  <CR> 

(2)  Choose  ".Set"  Irom  the  .Set/Adjust  '  menu  Select  with  <CR  >  (The  "Adjust" 
selection  is  used,  for  example,  when  the  charge  on  a  system  ol  biased  conductors 
is  to  be  initiali/.ed  I 

The  conductor  potential  screen  now  appears  as  shown  in  figure  2  4  The  "BC 
Type"  column  is  a  toggle  between  the  values  "I'ixed  Potential"  and  "Biased 
Potential"  If  you  wish  a  fixed  potential  value,  choose  "Fixed  Potential"  and  enter 
the  value  in  the  third  column.  If  you  wish  to  bias  the  coiuluclor  relative  to 


another,  choose  "Biased  Potential"  and  enter  the  relative  \alue  and  the  conductor 
number  to  which  it  is  biased.  (The  potential  dillcrcnce  of  biased  condueto  s  will 
remain  fixed  when  conductor  potentials  change.) 


Cond.  BC  Type  Value 

1  Fixed  Potential  O.tXXlE+OO 

2  Fixed  Potential  ().(XX)E+(K) 

Fixed  Potential  ().()(X)E+(X) 

4  Fixed  Potential  ().(XX)E+(Kt 


Biased  From  i 

( 

i 

0  I 

0  i 

! 

0  i 

_ I 


Figure  2.4.  Conductor  potential  screen 

(4)  After  achieving  the  correct  array  of  insulator  potentials,  signify  cmnpletion  \Mth 
the  <ESC>  key.  You  will  be  returned  to  the  main  IP.S  pull-down  menu. 

(5)  Choose  "Write"  to  enter  your  potential  conditions  mio  the  DynaPAC  DataBu'^e 
files. 

You  may  now  wish  to  change  some  insulator  boundary  conditions 

(,6)  Select  "Insulators"  from  the  pull-down  menu. 

(7)  You  are  given  a  screen  to  specify  the  range  of  conductor  numbers  and  the  range 
of  material  numbers  to  which  your  changes  will  apply.  (Only  insulating  materials 
will  be  affected.)  Alter  the  screen  to  suit  your  wishes,  and  signify  completion 
with  <ESC>. 

(8)  Select  "Fixed"  to  apply  fixed  potentials,  or  "Roat  "  to  apply  fixed  electric  field 
(Non-zero  fixed  electric  field  may  be  specified.)  The  remaining  choices  are  "Nor¬ 
mal"  (to  further  narrow  the  range  of  cells  to  which  the  condition  is  applied)  and 
"Edit"  (to  specify  condition  cell  by  cell  for  all  cells  in  the  range) 

(9)  Enter  the  value  of  potential  or  electric  field.  Signify  completion  with  <ESC>. 

(10)  To  see  the  result  of  your  entries,  select  ".Show  Table  "  The  table  displays  the 
matrix  of  conductor  numbers  and  material  names  with  the  iPS  '.;iecification  for 
each  matrix  cell.  The  entries  are: 


Table  Entry 

Meaning 

(Blank) 

i 

(Real  Number) 

FLOA 

MIX 

No  such  surfaces 

Fixed  potential  value  ; 

Fixed  electric  field  (of  some  value) 
Surfaces  have  differing  spe-cifications. 

(II)  Choose  "Write"  to  enter  your  potential  conditions  into  the  DynaPAC  DataBase 
files. 

IPS  writes  a  file,  "ipsout."  with  a  complete  cell  by  cell  specification  of  the  initial 
potentials. 
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2.4.2.  Using  DynaPre  1  Potent)  to  Create  Potential  Solver  Input 

Select  "Potent"  From  the  top  menu  of  (he  DynaPre  module.  Select  "Edit  Scnpt" 
from  the  resulting  pull-down  menu.  The  potential  solver  input  screen  shown  in  figure 
2.5  will  appear.  Modify  the  input  screen  to  suit  your  needs.  F'.xit  from  this  screen 
with  an  <ESC>.  Then  select  "Make  Script”  to  write  a  fully  commented  input  deck  for 
the  potential  solver.  Subsequent  modifications  to  the  input  deck  arc  conveniently  done 
with  the  text  editor. 

Here  is  a  discussion  of  the  potential  solver  input  parameters. 

Run  option: 

This  toggle  field  takes  the  values  "New"  or  "Continue."  "New"  must  be  specified 
for  the  first  potential  run  on  a  DynaPAC  file  set.  "Continue"  can  be  specified 
(and  is  usually  preferable)  on  any  subsequent  run,  even  if  IPS  has  been  used  to 
change  surface  boundary  conditions. 

.Algorithm: 

The  potential  solver  can  operate  in  "32  Node"  (continuous  electric  field)  mode  or 
"8  Node"  (trilincar  potential)  mode.  Other  modules  do  not  support  the  "8  Node" 
potentials,  although  the  potential  soliituin  is  much  faster.  It  has  not  been  demon 
strated  that  solving  potentials  with  the  8  Node  algorithm  and  then  switching  to 
the  32  Node  algorithm  saves  any  time 


^  Exit  CadDyn  IPS 


winterm 


Part Gen  Tracker  Help 


JJl 


Potential  Solver  Parameters 


Comment: 

1^ 

Run  option:  NEW 

Algorithm  32_NODE 

Problem  type;  LAPLACE 

Debye  Scaling; 
Solution  Mix; 
Moke  E  f  fee  t : 

LOCAL 

O.OOQE+00 

OFF 

Convergence  criteria: 

Environment: 

T imer : 

1 

Maxits  =  10 

Temp  =  I.OOOE-OI 

Diagnostics : 

RdRmin  =  l.OOOE-03 

Dens  =  1 . 000E+ 1 1 

Initial 

1 

RMSmin  =  l.OOOE-02 

Debye  =  7.434E-03 

Final : 

1 

Maxitc  =  40 

10  files  : 

SCG 

1 

PotCon  =  2.000E+00 

PS_input  :  ps_in 

Screen 

0 

DebLim  =  2.000E+00 

PS_output:  ps_out 

Special 

0 

Conv.  Effect:  ON 

Time  parameters: 

Inter  f 

1 

Selected  Grids; 

From  tt  0  to  t$  0 

Start  =  O.OOOE+00 

Rise  =  O.OOOE-hOO 

Fall  =  I.OOOE+30 

Wake 

1 

Figure  2.5.  Potential  solver  input  screen  from  the  DynaPre  preprocessor. 


Problem  Type: 

This  toggle  field  describes  the  space  charge  lunction  to  be  used  in  the  potential 

solver.  These  are  implemented  in  the  subroutine  rhwfp.t'.  Available  choices  are: 

LAPLACE  -  Zero  space  charge. 

LINEAR  -  Linear  (Debye)  screening 

NON  LINEAR  -  .Standard  equilibnum  space  charge  tormula. 

FROZEN  ION  -  (on  density  is  everywhere  equal  to  the  ambient  plasnui  density 
Electron  density  is  exp((|>/0)  for  negative  potentials. 

TRACK  ION  -  Ion  density  is  taken  from  particle  tracking  results  (array 
RHO  Ion).  Electron  density  is  exp(({»/(3)  for  negative  potentials. 

EXPLICIT  -  Both  electron  and  ion  densities  are  taken  frotn  particle  tracking 
results  (array  RHO  Ion).  ("Full  PIC.") 

SHEATH  WAKE  Ion  and  electron  densities  are  calculated  with  orbital  motion 
(i.e..  the  wake  of  the  spacecraft  sheath  and  its  effect  on  currents)  taken  into 
account. 

SPECIAL  -  Developmental  space  charge  option.  Presently,  this  option  is  similar 
to  the  SHFIATH  WAKE  option,  but  uses  tracked  electron  densities  (array 
RHO  Elcc)  m  positive  potential  regions. 

Convergence  criteria: 

Normally  only  the  two  parameters  "Maxits"  anil  "RMSmin"  are  varied  from  their 
default  values. 

Maxits  -  The  maximum  number  of  major,  or  "space  charge"  potential  iterations  to 
be  performed 

RdRmin  -  The  value  of  the  "RDOTR"  parameter  below  which  the  potential  will 
be  considered  converged. 

RMSmin  -  The  value  of  the  "RMSERR"  parameter  below  which  the  potential  will 
be  considered  converged. 

Maxitc  -  The  maximum  number  of  minor  iterations  within  each  conjugate  gra¬ 
dient  solution 

PotCon  -  The  number  of  orders  of  magnitude  that  the  RDOTR  drops  within  each 
major  iteration  before  it  is  considered  converged 

Debl.im  -  The  number  <if  Debye  screening  lengths  allowed  jvi  /one  Fhe  various 
space  charge  treatments  limit  the  amount  of  space  charge  in  a  /one  in  accoi 
dance  with  this  parameter. 

Conv.  Effect  Whether  the  analytic  trajectory  convergence  treatment  is  used  in 
the  NON  LINE-AR  problem  type. 


fi) 


Selected  (Jrids. 


It  is  possible  to  apply  the  potential  solver  to  a  subset  ol  the  DynaPAC  gnds 
This  option  has  not  been  fully  tested. 

Debye  Scaling 

The  "DebLim"  parameter  (see  above)  may  be  applied  based  on  the  local  grid 
spacing  or  on  the  primary  grid  spacing. 

Solution  Mix 

For  problems  which  do  not  converge  well  this  parameter  allows  admixture  of  the 
solution  from  the  previous  major  ("space  charge")  iteration  at  the  end  of  the 
current  major  iteration. 

Wake  Effect 

This  option  is  to  be  turned  on  to  take  account  of  spacecraft  motion.  It  brings 
down  a  sub-menu  to  allow  the  spacecraft  velocity  to  be  specified.  On  a  "NEW 
potential  run,  the  wake  of  the  (uncharged)  spacecraft  surfaces  will  be  calculated. 
("Neutral  Approximation.")  A  separate  module  is  used  to  calculate  the  wake  of 
the  spacecraft  sheath. 

Environment 

A  single  maxwellian  environment  is  specified  in  terms  of  plasma  temperature 
(eV).  plasma  density  (m  and  Debye  length  |rti|.  Only  two  of  these  parameters 
are  independent.  Debye  length  will  be  automatically  recomputed  following 
change  in  temperature  or  density.  Density  will  be  automatically  recomputed  fol¬ 
lowing  change  in  Debye  length. 

lO  Files 

The  PS  input  name  is  the  name  of  the  file  to  be  written  by  the  "Make  Script" 
selection  from  the  Potent  pull-down  menu.  PS  (lutput  presently  has  no  effect 

Time  Parameters 

For  time -dependent  problems  conductor  potentials  can  be  modified  (relative  to 
IPS  settings)  by  a  pulse-shape  factor.  This  is  done  in  subroutine  setend.F.  The 
expression  cunently  used  is 

P(t)  -  PCond*|  1  -exp(-(t-tl  )/t2)l*exp(-(t-t  1  )/t3) 
where  tl  is  start  time  Isecondsj 

t2  is  rise  time  (seconds] 
t.f  is  fall  time  |scconds| 
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(A  rise-time  of  zero  indicates  that  no  time-dependent  factor  is  to  be  applied.)  The 

time  parameter,  t.  is  advanced  by  the  Tracker  module. 

Timer  and  Diagnostics 

These  parameters  govern  optional  printtnit  from  various  portions  of  the  potential 

solver. 

2.4.3.  Operation  of  the  Potential  Solver 

In  this  section  we  discuss  the  operation  of  the  potential  solver  from  the  point  of 
view  of  monitoring  its  progress. 

It  is  not  recommended  that  the  potential  solver  be  run  from  within  DynaPre.  If 
Dynapac. setup"  has  been  sourced,  run  the  potential  solver  with 

Potent  <  ps  in  >  ps  out 

t)therwise,  use 

SDYN.APAC/Potent  eputype  <  ps  in  >  ps_out 
where  eputype  is  either  SUN’4  or  IRIX. 

A  potential  solver  run  consists  of  an  imtiali/.ation  pha.se,  a  number  of  major 
("space  charge")  iterations  of  the  potential  solution,  and  a  brief  exit  phase.  In  the  ini 
tializ,ation  phase,  the  input  parameters  are  read  and  echoed,  and  the  DataBase  informa¬ 
tion  is  processed.  Grid  information  (mesh  size  and  sheath  potential)  is  printed.  The 
grid  interface  pairs  list  is  formed  and  printed.  Conductor  potentials  are  printed. 

At  the  beginning  of  each  major  potential  iteration  the  space  charge  function  and 
its  derivative  are  evaluated  cell  by  cell,  and  the  conjugate  gradient  process  is  initial¬ 
ized.  The  "initial  rdotr"  is  a  measure  of  the  current  error  in  the  potential  solution.  For 
most  cases  the  "initial  rdotr"  should  decrease  monotonically  beyond  the  first  few  majoi 
iterations. 

During  the  conjugate  gradient  process  (governed  by  subroutine  psscg.F)  the 
"rdotr"  parameter  is  printed  lor  each  minor  iteration.  This  parameter  should  generally 
decrease,  but  for  most  cases  will  not  decrease  monotonically,  Ihe  command  "grep 
rdotr  ps  out"  is  a  good  way  to  check  on  the  progress  of  the  potential  solver.  The  con- 
lugate  gradient  process  is  deemed  converged  when  the  criteria  discussed  above  (usually 
twd  orders  of  magnitude  decrease  in  "rdotr")  arc  satisfied 

We  are  now  ready  to  conclude  the  major  iteration.  The  most  i  nie-consuming 
task  is  to  calculate  and  update  surface  electric  liekls,  which  are  held  fixed  during  the 
conjugate  gradient  process  This  requires  a  full  matrix  operation  (performed  by  sub¬ 
routine  ecoprd.F),  as  the  electric  field  value  is  related  to  the  residual  for  the 
corresponding  potential  The  held  updates  (subject  to  several  rest rict ions)  are  per- 
lormeil  bv  subroutine  eupd.it  F'  The  new  surface  potentials  and  fields  arc  printed  The 


ditTerencc  between  the  new  and  previous  potentiai  solulioiis  are  expressed  a  ruoi 
mean-square  errors,  and  are  printed  grid  by  grid  (  "grep  RMS  ps  out  )  and  oxeral! 
("grep  rniseiT  ps  out").  If  the  root-rnean-square  errors  remain  constant,  solution- 
mixing  may  ameliorate  or  .solve  the  problem.  (A  particularly  large  "RMS  Error"  tor  a 
particular  grid  may  or  may  not  indicate  a  problem  in  obtaining  a  solution.) 

The  Potential  .Solver  is  ready  to  conclude  when  the  requested  nurnlx'r  of  ma|or 
iterations  have  been  performed,  or  when  the  "rmserr'  has  been  reduced  below  its 
minimum  value.  By  default,  the  central  Z-slice  potentials  and  electric  tields  are 
printed  for  each  grid.  In  these  printouts  three  values  are  printed  for  each  grid  point; 
the  potential  value  appears  at  the  print  position  corresponding  to  the  grid  point;  the  x- 
direction  potential  gradient  [volts  per  melcrj  is  printed  to  the  right  of  the  grid  point, 
and  the  y-direction  potential  gradient  above  the  grid  point.  The  /-direction  potential 
gradient  does  not  appear. 

2.4.4.  Using  Scanner  to  Display  Potentials 

The  Scanner  module  is  a  general  interactive  utility  to  extract  and  display  surface¬ 
cell  or  grid  structured  information  from  the  DataBase,  Its  main  use.  however,  is  to 
plot  potentials. 

The  top  menu  biu"  ol  Scanner  contains  the  selections  "E-xit."  "Print."  and  Plot 
With  the  "Print"  selection,  virtually  any  surface  or  spatial  array  known  to  the  DataBase 
Manager  can  be  printed.  The  print  facility  is  primarily  used  for  diagnostic  purposes, 
and  is  fairly  self-explanatory. 

When  "Plot"  is  selected,  a  pull-down  menu  is  selected  with  the  options  '  Exit.' . 
"Data  Info,".  "Edit  Plot,"  "Make  Plot",  "Show  Plot."  Each  of  these  options  (excluding 
"Exit")  should  be  selected  at  least  once  prior  to  .selecting  any  of  those  below.  The 
options  chosen  for  "f)ata  Info"  and  "Edit  Plot"  are  saved  (upon  exiting  .Scanner)  in  the 
tile  "prefix. scan  "  If  you  are  making  a  few  different  plots,  you  may  wish  to  make  a 
copy  of  this  file  for  each  plot,  and  copy  it  back  to  "prefix  scan"  to  efficiently  regen¬ 
erate  your  plot  options. 

The  "Data  Info"  selection  allows  you  to  choose  the  quantity  to  be  plotted  and  the 
norma!  direction  to  the  plot  plane  The  "Data  Name"  field  is  a  toggle  field,  currently 
having  (he  options  POT  Grid  (potentials).  RflO  Gl  (ion  densities  from  wake  shadow¬ 
ing  calculation).  RHO  Elec  (tracked  electron  densities),  RHO  Ion  (tracked  ion  densi 
ties),  NONE  (plot  gridding  only),  and  OTHER  (allowing  user  specification  of  another 
name). 

The  T'klil  Plot"  selection  gives  the  user  a  great  deal  of  control  over  the  plot. 
Most  of  the  options  are  reasonably  self-explanatory.  Be  sure  that  the  "Window 
Manager"  and  "Plotter"  toggle  fields  are  correctly  set.  Selections  under  the  "Options" 
title  each  bring  up  a  sub-menu;  "GRID  I.IMITS"  allows  the  user  to  specify  the  region 
of  space  to  be  plotted  and  the  units  for  the  axis  labels  "GONTOfiR  f.EiVELS"  allows 


spccilication  of  up  to  24  contour  levels  (or  colt)r  boundaries)  "CONTOUK  MARKS" 
provides  for  labeling  of  contour  lines.  "LABiil.S  Ac  INt’l.UDES"  controls  the  labeling 
and  appearance  of  the  plot.  "OTHER  Ol’TIONS”  at  present  contains  obsolete  or 
little-used  choices. 

The  "Make  Plot"  selection  generates  the  plot  and  writes  it  on  the  disk  file 
"fort. 2."  If  there  are  unvieued  plots,  the  new  plot  will  be  appended  to  the  "fort  2”  (ilc. 
If  there  are  no  unvievved  plots,  any  existing  "fort  2"  will  be  overwritten 

The  "Show  Plot"  selection  brings  up  a  shell  to  run  the  chosen  plot  display  inter 
face.  The  tile  "fort  2"  will  remain  after  exit  from  the  .Scanner  module,  and  can  be 
displayed  or  printed  using  any  plot  interface  program  in  stand  alone  mode 
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2.5.  (k'ncratin^  and  Trackins  Particles 

Macropeulicks  may  be  used  in  DynaPAC  lor  a  number  ol  purposes  includin': 

(1)  Studying  and/or  displaying  representative  particle  lra|eciones; 

(2)  Calculating  surface  currents  arising  from  sheath  currents, 

(3)  Studying  wake  structure; 

(4)  Calculating  steaily  state,  sell  consistent  charge  ticnsities; 

(5)  Calculating  time-dependent  charge  densities  and  surface  currents 

To  maximi/e  Ilexibility,  the  particle  generator  and  the  particle  tracker  have  been 
developed  as  two  separate  modules  fhe  particle  generator  defines  parti*  le  species  and 
generates  particles  throughout  volume,  aUmg  a  sheath  surface,  along  a  con  line,  in 
accordance  with  external  input,  or  along  magnetic  field  lines.  The  particle  tracker  will 
compute  the  motion  of  all  or  a  subset  of  the  particles  for  a  maximum  time  or  end 
cycling,  recording  surtace  currents,  and  (optionally)  plotline  irajcctories.  accumulating 
steady-state  charge  density,  or  calculating  new  charge  density  at  the  updated  time 
After  the  panicle  tracker  is  run.  the  particle  files  are  left  wi:h  updated  particle  posi¬ 
tions  and  velocities,  and  the  time  and  cycle  number  is  updated  in  the  DataBase. 

Note  that  there  are  obvious  relationships  between  the  Potential  .Solver  module's 
"Problem  Type"  space  charge  options,  the  Panicle  Generator  module’s  "Panicle  Tvfx;" 
options  and  the  Panicle  Tracker  module  s  "Process”  options  It  is  prei  urc  to  lay 
down  firm  rules  about  which  options  correspond,  but  the  user  should  pay  attention  lo 
the  physical  processes  modeled  and  the  data  anays  generated  or  required  in  planning  a 
calculation  strategy. 

2.5.1.  U.sing  Dynal’re  (Fas’.Jenl  fo  Create  Particle  (Jenerator  Input 

The  particle  generator  input  screen  is  obtained  with  the  'PartGen/Edit  Senpt" 
selection  from  DynaPre.  ft  is  used  to  obtain  a  fully  commented  input  deck  for  the 
PartGen  module.  I  he  text  editor  may  be  used  subsequently  to  modify  the  input  deck 

The  main  screv..  ha.s  .selections  for  the  magnetic  field,  the  number  of  species,  and 
the  name,  charge,  and  mass  of  each  species.  IXtnsity  and  temperature  will  be’  the 
values  used  on  the  most  recent  potential  solver  run. 

Further  description  of  the  particle  generation  process  depends  on  the  '  Particle 
Type"  toggle  field.  The  choices  are; 

DKFADLT 

(Jenerate  particles  (S  per  /one)  representing  a  uniform  distribution  of  charge 


SHKATH 


Generate  particles  representing  sheath  ciirienis.  A  suhscicen  will  appear  to 
enter  the  sheath  potential 

CONTOUR 

Generate  particles  along  a  contour  line.  (This  is  usually  clone  in  order  to 
display  particle  trajectories.)  A  subscreen  will  appear  to  enter  the  contour 
pcMential  and  the  normal  and  ot'lset  of  the  cut  plane  (in  pnmary  grid  units). 

EXTERNAL 

Read  initial  particle  data  from  an  external  tile.  A  subscrecn  appears  for  the 
file  name.  I  he  external  file  is  read  by  subroutine  redext f  The  file  formal 
is 

POSITION  \  y  / 

niRECTION  dx  dy  d/ 

ENERGY  energy 

where  the  position  is  in  primary  grid  units,  the  direction  need  not  be  normal- 
i/.ed.  and  the  energy  is  kinetic  energy  in  electron  volts  Alternatively,  the 
keyword  E  TOTAl.  may  be  used  to  give  the  total  panicle  energy.  Each 
such  three  card  sequence  results  in  definition  of  a  pculicle.  A  unit  particle 
weight  is  assigned,  on  the  assumption  that  this  input  is  used  only  for  trajec¬ 
tory  plotting  purposes. 

B  FIELD 

The  B  FIELD  option  is  used  to  generate  particles  where  magnetic  field  lines 
enter  the  computational  space.  The  subscreen  is  used  to  enter  the  particle 
energy,  the  (square  root  of  Ihe)  number  of  particles  Ic  r  each  boundary  sur 
face,  and  the  magnetic  field. 

2.5.2.  U.sing  DynaPre  ri  rackerl  to  Create  Particle  Tracker  Input 

The  particle  tracker  input  screen  is  obtained  with  the  '  Tracker/Edit  Script"  selec¬ 
tion  from  DynaPre.  It  is  used  to  obtain  a  fjlly  commented  input  deck  for  the  Tracker 
module  The  text  editor  may  he  used  subsequently  to  modify  the  input  deck 

The  main  screen  has  selections  for  Ihe  magnetic  field,  the  number  and  selection  of 
species  to  be  tracked,  the  maximum  trajectory  lime,  other  limits  concerning  tracking 
fidelity,  and  limits  associated  with  particle  origins  and  trajectory  display  There  is  a 
slot  reserved  for  the  particle  pushing  algorithm,  though  at  present  only  one  algorithm 
Is  implemented  (in  subroutine  movpar.FV 


The  use  1(1  which  particle  trajectory  intorriiatioii  is  put  (usually  in  subroutine 
savinf.F)  is  governed  by  the  "Proeess"  parameter  (which  is  a  toggle  tield)  (.  uriem 
choices  for  the  "Process"  parameter  are: 

SPACE  CIIAR(;E 

This  is  the  "PIC"  option.  After  each  particle  has  been  tracked  for  the 
requested  time,  its  charge  is  shared  to  the  nodes  of  its  current  grid  cell  The 
charge  is  stored  in  the  KIK)  Ion  array  Special  cells  and  surface  elements 
are  treated  correctly. 

TRAJEC  rORIES 

This  is  the  option  to  track  particles  for  the  purpose  of  either  displaying  their 
trajectories  or  calculating  surface  currents.  A  subscreen  comes  up  to  select 
the  direction  normal  to  the  plot.  (Only  one  plot  is  produced.)  To  obtain  a 
plot,  the  "Plot  Limit"  option  (see  below)  must  be  specified;  otherwise  plot¬ 
ting  will  be  turned  off. 

PI.OT  PARIICEES 

With  this  option  no  tracking  is  performed,  but  the  current  positions  of  the 
particle  are  plotted.  A  subscreen  comes  up  to  select  the  direction  normal  to 
the  plot.  (Only  one  plot  is  produced.) 

TKI  CHARfJE 

This  is  the  option  to  accumulate  space  charge  over  the  entire  trajectory 
length.  The  space  charge  is  not  shared  to  nodes,  but  is  accumulated  as  an 
element  centered  array,  KHO  Elec.  Charge  accumulates  in  empty  cells  only 
Special  cells  and  surface  elements  are  not  treated  correctly. 

Other  particle  tracker  options  include: 

Tracking  Time  -  The  maximum  time  (seconds)  a  particle  is  to  be  tracked.  For  time 
dependent  problems,  the  time  variable  will  be  incremented  by  this  amount 

Max  Dx  -  The  maximum  distance  (in  local  grid  units)  a  particle  moves  during  a  sub¬ 
step.  Substep  distances  may  be  far  less  than  this  value.  Electrons  gyrating  a  a 
magnetic  field  move  only  for  a  fraction  of  the  cyclotron  period.  Slow  ions  may 
move  less  than  "Max  Dx"  in  the  tracking  time 

Max  .Steps  -  The  maximum  number  of  substeps  per  particle  per  iteration. 

Max  iterations  -  The  maximum  number  of  cycles  through  the  list  of  grids. 

Saving  Interval  -  How  often  potentials  arc  saved  (by  subroutine  spaini.F)  in  time 
dependent  problems. 


Particle  Limit  -  Limits  (in  primary  grid  units)  t)(  initial  particle  locations  Particles  on 
ginating  outside  these  limits  are  ignored  II  not  set  (all  zeroes)  all  pani«  les  uil! 
be  tracked. 

Plot  Limit  -  Limits  within  which  particles  or  trajectories  are  pioiied  It  not  set  (all 
zeroes)  no  plot  will  be  produced  by  the  I  K.AJPX'TOK  Y  prtKCss 

Number  of  Species  -  Number  of  species  to  be  tracked 

Species  Code  -  Which  s|Kcies  to  be  tracked  (corresponding  to  species  definitions  in 
particle  generator  input) 

B  field  -  Magnetic  field. 

Mixing  Factor  -  For  the  TRJ  CHARGF  process,  the  fraction  of  the  previous 
RHO  Elec  values  to  mix  with  the  newly  computed  values 

TR  input  -  Name  of  file  produced  by  "Make  Script"  selection. 

TR  output  -  (No  effect) 

Tracker  Diag  -  Diagnostic  level  for  Tracker  subroutines. 

Dynalib  Diag  -  Diagnostic  level  for  Dynalib  subroutines. 

Timer  Level  -  Frequency  of  epu  time  monitoring 

2.5.3.  Operation  of  (he  Particle  (Jenerator 

In  this  section  wc  briefly  describe  how  the  particle  generator  operates  for  each 
"particle  type,"  pointing  to  some  of  the  key  subroutines  in  the  process.  In  general, 
subroutine  inivel.F  is  used  to  assign  initial  velocities  to  particles.  Subroutine  convrt.F 
converts  volume  weights  to  charge  or  current  weights. 

For  the  "Default"  panicle  type  (IFType=l,  uniform  density  throughout  space), 
particle  generation  is  under  the  control  of  subroutine  genpar.F.  Actual  particle  genera¬ 
tion  is  done  by  subroutine  fndpar.F.  Particle  weights  are  proportional  to  the  shadow¬ 
ing  factors,  RHO  Cl  Particle  weights  are  charge  divided  by  the  permittivity  of  free 
space,  or  [volts  m  "[, 

Fur  the  "SHEATH"  particle  type  (lPType=2,  particles  generated  on  potential  con 
tour  surfaces)  particle  generation  is  under  the  control  of  subroutine  genpa2.F.  Actual 
particle  generation  is  done  by  subroutine  sihpar.F.  SthPar  varies  its  procedure  if  mag 
netic  effects  dominate,  as  determined  by  comparing  the  Larmor  radius  to  the  outer  grid 
spacing.  Subroutine  do/one  F  determines  the  panicle  positions  and  .irea  weights,  also 
eliminating  high  field  or  bipolar  zones.  For  non-magnetic  cases,  subroutines  inivel.F 
and  dflux.F  detennine  panicle  velocity  and  current.  For  magnetic  cases  subroutine 
niagsth  performs  this  function,  assigning  a  "cos  0"  factor  to  account  for  particles  trav¬ 
eling  along  field  lines  to  reach  the  sheath  surface  Panicle  weights  are  amperes. 

For  the  "CONTOUR"  particle  type  (IPType=-^).  particle  generation  is  under  the 
control  of  subroutine  genpa.^  F  Actual  particle  generation  is  done  by  subroutine 
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cntpar.F*. 

For  the  "F\ri:RNAL"  purtieic  lyjK  tlPType=4).  |)ailiele  generation  is  under  the 
eoiitrol  ot‘ subroutine  genpa4.F.  Aetual  pailiele  generatit>n  (l  e  ,  reading  the  input  tile) 
is  done  by  subroutine  redext.F. 

For  the  "B  FIELD  '  particle  type  (IPType=5).  particle  generation  is  under  the  con 
trol  of  subroutine  genpa5.F  Particles  are  generated  only  on  the  surface  of  the  outer¬ 
most  grid.  Actual  particle  generation  is  ilone  by  subroutine  bparts.f  Particle  weights 
are  amperes. 

2.5.4.  Operation  of  the  Particle  Tracker 

The  particle  tracker  operates  on  the  prefix. PTn  tiles,  which  are  organized  in  pages 
of  lUOO  panicles  each.  Track  is  kept  ol  how-  many  particles  on  each  page  belong  to 
each  grid.  An  "iteration"  consists  of  Itioping  through  each  grid,  and  within  that 
through  the  panicles  belonging  to  that  grid.  .Subroutine  trkpar.F  supervises  the  track¬ 
ing  of  each  particle  by  subroutine  movpar.F.  Each  particle  is  tracked  until  one  of  the 
following  conditions  occurs: 

(1)  The  particle  strikes  the  spacecraft; 

(2)  The  particle  is  found  in  a  grid  other  than  the  current  grid  or  its  parent,  in  which 
case  it  is  transferred  to  the  new  grid. 

The  particle  exits  the  computational  space; 

(4)  The  trajectory  time  reaches  the  requested  particle  tracking  time; 

(.5)  The  number  of  substeps  exceeds  the  maximum  substep  number; 

(6)  For  the  TKJ  CHARGE  process,  it  the  particle  is  loiind  outside  the  "^heath."  (i  e  . 
in  a  region  where  it  has  very  low  kinetic  energy). 

Trajectories  stopped  for  conditions  (2)  or  {5)  above  are  picked  up  during  the  next  itera¬ 
tion. 

After  processing  each  grid  Tracker  prints  the  particle  number  and  total  weight  lor 
each  of  several  particle  categories: 

(1)  "new  particles"  is  the  set  of  particles  with  which  Tracker  started; 

(2)  "partially  tracked"  refers  to  particles  which  have  not  hit  the  spacecraft  or  left  the 
primary  grid.  These  particles  may  or  may  not  yet  have  been  tracked  for  the 
requested  particle  tracking  time. 

(.'^)  "dead"  particles  are  those  which  struck  the  spacecraft 

(4)  "went  off  primary  grid"  counts  those  particles  which  have  exited  (he  computa¬ 
tional  space. 

(.5)  "trapped"  is  presently  a  null  category,  as  there  is  no  criterion  for  (rapped  particles 
Trapped  particles  appear  as  partittliy  tracked  " 
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(6)  "unknown  status'  is  another  category  of  particles  with  no  known  means  of 
identification. 

When  Tracker  is  finished  with  the  requested  number  of  iterations,  it  pnnfs  the 
total  current  to  the  spacecraft  and  a  table  appoilioning  ihai  current  by  material 
and  conductor.  Also,  a  "Hit"  tile  is  printed  listing  each  particle  that  struck  the 
spacecraft  with  its  vital  statistics  which  may  be  used  for  further  processing.  The 
print  and  format  sialemeiits  for  this  file  may  fv  found  in  subroutine  wripar  F 
If  the  process  invoKcs  traiectory  plotting.  .i  "foil  2  "  file  will  K-  created  containing 
the  plot. 
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2.6.  Using  the  DynaPost  Utilities 

DynaFost  is  a  u.Ncr  interface  to  the  ('urrenis  and  ObiFoil  postprocessors 

The  Currents  postprocessor  is  intended  priinanly  lor  iinic-dependciit  nrohlems 
Hy  reading  data  from  the  prefix. CT IK  life,  it  can  reproduce  the  table  of  surface  em  n  ni 
(sorted  by  conductor  and  material  number)  printed  at  the  conclusion  ol  the  Irackct 
module.  It  can  also  create  plots  of  one  or  more  of  the  table  entries  as  functions  of 
time. 

The  ObjPotl  postprocessor  allows  the  user  to  display  the  spacecraft  and  the  results 
for  spacecraft  surfaces  Plots  produced  are 

(1)  Object  outline  t  "wireframe” )  plot.  ct»lor  coded  by  niatcnal 

(2)  Object  with  surface  potential  contour  lines  drawn 
(.^)  Object  with  surface  cells  shaded  by  matenal 

(4)  Color-coded  object  surface  potentials 

(5)  Color-coded  surface  normal  electric  fields. 

(6)  Object  with  surface  cells  shaded  by  conductor  number, 

(7)  Object  with  surface  cells  shaded  by  incident  current 

(8)  Color-coded  incident  Hux  densities 

When  ObjPotl  is  chosen  from  the  top-level  DynaPost  menu,  the  user  is  presented 
with  a  screen  from  which  he  can  vary  input  options  to  ObjPotl.  create  an  input  file,  run 
(yes.  this  one  actually  works!!)  ObjPotl.  and  view  the  resulting  plot  file  General 
options  include  a  title  for  the  plot,  a  vector  tunnormali/ed)  f  rom  the  spacecraft  toward 
the  viewer,  and  the  plot  display  interlace  to  be  used.  By  default,  surface  cell  numbers 
appear  on  plots  (I).  (.1),  and  lb)  above;  toggles  to  "Aid."  or  'NONE:'"  control  the 
appearance  of  these  numbers.  E'or  each  plot  a  "YK.S/NO"  toggle  is  a\ailable  to  deter¬ 
mine  whether  the  plot  is  made,  in  addition  plots  with  no  information  (e  g.,  uniform 
potential)  will  not  be  made.  For  the  potential,  field,  current  and  flux  plots,  the  range 
of  the  color  table  may  be  restricted;  by  default  the  range  with  be  the  range  of  actual 
values  found. 

The  run  controls  appear  in  the  lower  left  of  the  screen  "MAKE  SCRIPT"  writes 
a  fully  commented  input  deck  as  "'objpotl.in.""  "RUN  .SC  RIPT"  issues  (in  either  fore¬ 
ground  or  background)  the  command  "'ObjPotl  <  objpotl  in  >  objpotl.out."  creating  a 
"fob. 2"  plot  file.  ""SHOW  PLOT"  spawns  a  shell  to  run  the  selected  plot  display  inter¬ 
face. 

On  exit  a  plot  tile  "  fort. 2”  will  remain  for  disposition  by  any  of  the  plot  display 
interfaces,  and  the  input  and  output  files  from  ObjPotl  will  remain.  If  the  incident  flux 
plot  has  been  made,  the  area,  current,  and  flux  to  each  surface  cell  will  be  listed  in  the 
tile  "Objpotl Currents,"" 
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2.7.  Time-Dependent  Problems 

DynaPAC  is  intended  for  time-dependent  prc^blems  However,  it  must  be  con 
fessed  that  only  one  complex  time-dependent  problem  (involving  mobile  ions  and 
bru'ometric  electrons)  has  ever  been  run  It  k  to  he  expected  that  any  application  not 
entirely  similar  to  that  one  will  require  some  code  nuxiitication  In  this  section  we 
remind  the  user  o(  the  time-  dependent  leatuics  ot  DynaPAC.  niosi  of  which  were 
mentioned  above  Remember  that  the  bulk  ol  a  time  dependent  problem  consists  ot 
alternate  Potent-Tracker  runs. 

2.7.1.  Setting  Initial  Potentials 

In  DynaPre/IPS,  potentials  may  be  set  to  fixed  values,  or  they  may  be  set  to 
"Bias'  values  and  the  "Adjust'  options  used  In  the  f'ortner  case  ( lived  values),  the 
pulse-shape  options  ot‘  Potent  may  be  used  to  create  a  time  depeinlent  multiplier  lor 
the  fixed  potentials.  In  the  latter  case  ("Adjust"  options)  the  overall  object  potential 
may  change  due  to  currents  incident  on  the  spacecraft 

2.7.2.  Running  Potent 

Potent  should  be  run  with  appropnate  "Ihoblem  Type",  such  as  "IRACK  ION" 
or  HXPLICIT  "  However.  "f'ROZRN  ION"  or  "I.APl.AC'l: "  may  be  more  appropriate 
to  the  initiuli/ing  run  Remember  to  set  pulse  shape  parameters  if  ilesired.  and.  after 
the  first  run.  to  set  the  "New  or  Continue"  parameter  to  "Continue."  Also,  for  the 
"Continue"  runs  you  probably  will  want  to  set  "Maxits"  to  a  low  value  (such  as  two  or 
three). 

2.7..^.  Running  PartOen 

Normally,  for  time-dependent  problems  PartGen  will  he  run  after  the  initial  Potent 
with  the  "DKFAUIT"  panicle  type,  corresponding  to  a  uniform  charge  distribution 
There  is  currently  no  provision  for  refreshing  the  panicles  as  (hey  are  drained  from  the 
computational  space.  There  also  are  no  current  provisions  for  secmidary  particles  at 
surfaces  or  for  ioni/.ation  products. 

2.7.4.  Running  Tracker 

Tracker  is  the  key  module  for  lime  ilcpendeni  problems  With  the 
"SPACE  CHARGE"  process,  a  Traekcr  run  will  update  the  time  and  timestep  number, 
save  the  surface  currents  (m  the  prefix.CUR  file),  and  (in  accordance  with  Tracker 
input)  save  the  spatial  and  surface  potentials. 
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2.7.5.  Displaying  Results 

Scanner  can  display  any  of  the  poicniial  arrays  saved  by  ihe  Tracker  module. 
DynaPost  can  form  tables  and  plots  of  surface  currents  or  their  time  histories 

Tracker  may  be  used  with  the  "Pt.OT  PARTICLES  "  process  to  display  current 
particle  positions.  {This  is  usually  done  immediately  following  the 
Tracker/SPACE  CHARGE  run.) 


3.  DynaPAC  Developer’s  Manual 

This  section  is  intended  primarily  for  users  who  would  like  to  modify  or  extend 
the  capabilities  of  DynaPAC.  There  is  not  much  in  this  initial  version,  as  the  contract 
has  run  out  of  steam. 

3.1.  DynaPAC  Software  Structure 

DynaPAC  has  three  main  portions  contained  in  three  directories; 

The  "dynapac"  directory  contains  all  of  the  physics  modules,  as  well  as  many  of 
the  utilities.  This  directory  cosTesponds  to  the  DYNAPAC  environment  variable. 

The  "ScreenPkg"  directory  contains  the  low-level  routines  for  handling  the 
interactive  screens. 

The  "Plotters"  directory  contains  the  plot  display  interface  routines. 

In  this  manual  we  will  only  be  concerned  with  the  "dynapac"  directory. 

3.1.1.  DynaPAC  Directory  Free 

In  this  section  we  describe  contents  (d  the  SDYNAPAC  directory  and  its  sub¬ 
directories.  We  use  an  appended  slash  (!)  to  indicate  a  directory;  absence  of  a  slash 
indicates  a  tile. 

DynaPAC.. setup  -  This  is  the  file  establishing  environment  variables  and  aliases  for  a 
DynaPAC  user  or  developer. 

Makefile  -  This  should  tnake  all  the  DynaPAC  modules. 

bin/  -  Directory  for  DynaPAC  executable  modules  and  shell  scripts. 

dynamake  -  Shell  script  which  runs  DynaPAC  Makefiles. 

inputs/  -  Contains  screen  definitions  and  defaults  for  the  DynaPAC  interactive 
modules. 

lib/  -  Contains  DynaPAC  licraries. 

sre/  -  Contains  the  source  code,  object  files,  and  include  files  for  the  DynaPAC 
modules  and  libraries 

sre/AddGrid/  -  Code  to  add  an  outer  grid  to  a  DynaPAC  grid  structure, 
sre/DynaPost/  -  Code  for  the  DynaPost  module. 

Ntc/l)ynaPrc/  -  Code  for  the  DynaPre  module. 
sre/CiridTool/  -  Code  lor  the  GridTool  module 

sre/NisDyn/  -  Code  for  object  definition  interface  to  EMRC  DISPLAY-11. 
src/ObjPotl/  -  Code  for  the  ObjPotl  module. 
sre/PartGen/  -  Code  for  the  PartGen  module. 
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src/PatDyn/  -  Code  lor  the  PatDyn  (lUKiule. 
sre/PolDyn/  -  Code  tor  the  PolDyn  module. 
sre/Potent/  -  Code  for  the  Potent  module. 
src/Scanner/  -  Code  for  the  Scanner  module. 
src/Tracker/  -  Code  for  the  Tracker  module. 

src/cadlib/  -  Code  for  the  library  of  object  definition  interface  routines. 

src/currents/  -  Code  for  the  library  of  routines  used  in  the  DynaPosi/Currents  module. 

sre/dblib/  -  Code  for  the  DataBa.se  Manager  library  routines. 

src/dynalib/  -  Code  for  utility  routines  used  by  PartGen  and  Tracker. 

sre/ipslib/  -  Code  for  routines  used  by  DynaPre/IPS. 

src/.s3utils/  -  Code  for  the  low-level  utility  routines  used  by  all  DynaPAC  modulcN 
src/*/Dbx/  -  Source  code  resulting  from  processing  by  the  '  epp”  preprocessor. 

3.1.2.  DynaPAC  Software  Conventions 

There  are  many  DynaPAC  software  conventions.  Sometimes  they  are  followed. 
Other  times  they  are  not. 

Always  use  "syopen"  to  obtain  a  logical  unit  number  and  open  a  file.  Get  rid  of 
it  with  "syclos". 

3.1  J.  Modifying  DynaPAC  Subroutines 

Always  type  "../../dynamake"  to  remake  a  module. 

In  any  of  the  $DYNAPAC/src/mod«/eH«me  directories  you  will  find  source  files 
with  suffices  f,  .F.  and  .h. 

If  you  modify  a  .f  subroutine,  you  may  compile  it  in  the  usual  way,  then  type 
"../../dynamake"  to  remake  the  module.  If  you  do  not  fuanually  compile  it,  it  should 
be  compiled  by  the  Makefile. 

If  you  modify  a  .F  subroutine,  the  Makefile  should  run  the  "CPP"  processor  to 
produce  Dbx/name.f,  compile  Dbx/name.f,  and  remake  the  module. 

If  you  alter  a  .h  include  file,  the  Makefile  .should  reprocess  all  the  necessary  rou¬ 
tines.  Note  that  occasionally  include  files  are  shared  between  modules. 

If  you  add  a  routine,  be  sure  to  add  it  to  the  relevant  list  m  the  Makefile 

If  you  alter  a  library  routine,  you  may  manually  compile  it  and  add  it  to  the 
library,  or  you  may  type  "../../dynamake”  to  take  the  library  apart,  compile  the  routine, 
and  reconstruct  the  library. 


3.2.  Using  the  DataBase  Manager 

The  following  description  of  the  DataBase  Manager  appeared  in  the  interim 
Report,  SSS-DPR-90- 1 1973,  30  September  1990.  Check  in  the  17ynaFAC  modules  for 
further  examples  on  the  use  of  the  DataBase  Manager. 


3.2.1.  Data  Base  Library 

The  data  base  iibrary  provides  the  analysis  and  utilities  functions  necessary  to 
define,  manage,  store,  and  retrieve  data.  In  addition  to  the  three  callable  routines 
(dbinfoO,  dbdataO,  and  dbfileO)  normally  used  to  integrate  the  data  base  into  an  appli¬ 
cation,  a  stand-alone  driver  (dbtool)  is  also  available  to  interact  directly  with  the  data 
definitions,  data,  and  their  files. 

When  a  DyiiaPAC  module  is  linked  to  the  data  base  library,  the  following  rou¬ 
tines  provide  memory,  data,  and  file  management.  The  dbinfoO  routine  is  used  to 
define  data  base  building  blocks.  Typical  information  commands  define  grids,  grid 
nesting  structures,  data,  lists,  and  elements.  The  dbdataO  routine  is  used  to 
allocate/release  reusable  memory,  read/write  from  disk  files  to  memory,  and  initialize 
data  values  in  memory.  The  dbfileO  routine  provides  the  means  to  manage  the  data 
files  and  their  contents. 

The  input  to  these  routines  is  in  the  form  of  strings  containing  keyword  com¬ 
mands.  The  output  is  returned  to  the  calling  routine  via  common  blocks  within  spe¬ 
cialized  include  files.  The  calling  routine  saves  the  information  it  requires  in  its  own 
data  structures. 

The  data  base  can  be  thought  of  as  a  utility  driven  by  its  own  keyword  input 
language.  Although  some  defaults  are  built  into  the  data  base,  in  most  cases,  the  data 
required  to  perform  a  calculation  must  be  constructed  before  it  can  be  used  The 
dbinfoO  routine  is  used  to  define  data  and  the  dbfileO  routine  is  to  define  files. 

The  reader  is  referred  to  the  data  base  header  file,  DBhead,  for  a  complete  key¬ 
word  definition  and  current  implementation  status.  Some  examples  of  the  use  of  the 
library  routines  are  presented  below.  The  first  three  demonstrate  some  data  base  func¬ 
tions.  The  final  two  examples  show  how  a  calculation  using  data  in  the  data  base 
memory  storage  area  might  be  used. 
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3.2.2.  Example:  File  Detinition 


The  Task 

A  calculation  wants  to  open  a  set  of  DynaPAC  hies  that  have  the  prefix 
"Demo  Files"  and  then  close  them  when  it  finishes. 


The  Solution 

To  open  the  files,  use  the  command 

cal!  dbfiieC’open  prefix-Demo  Files") 

This  command  assumes  the  default  file  type  is  "DYNAPAC.  A  more  explicit  com¬ 
mand  would  be 

call  dbfiieC’open  prefix=Demo_Files  filejype=DyNAPAC") 

To  close  just  the  Demo  Files  file  set,  use 

call  dbfileC’close  prefix=Demo_Files") 

If  desired,  errors  opening  or  closing  the  files  can  be  checked  by  examining  the  error 
flags  in  tlie  dbfileO  output  include  file,  "odbfile.h”.  The  logical  flag,  Idferr  ,  is  TRUE 
if  ai'.  error  was  detected.  The  character  string,  "cdferr",  will  contain  an  error  message 
after  error  detection. 


3.2.3.  Example:  Definition  of  List  Dependent  Data 


The  Task 

Define  a  Surface_List  list  type  data  item,  "Surface  Cemroids ", 

The  Solution 

To  define  a  data  item  that  is  a  list,  both  the  list  type  and  the  data  item  must  be 
defined. 

Two  definitions  are  helpful  since  the  list  type  definition  depends  on  the  particular 
object  being  defined,  while  this  data  item  does  not.  Since  the  number  of  surfaces,  or 
length  of  Surfuce_List  is  not  given,  let’s  define  the  data  item  in  steps  first. 

Define  the  data  type  of  Surface_Centroids. 

call  dbinfoC’Define  Data=Surface_Centroids  Data_Type=LlST"} 

And  it  is  which  kind  of  LIST? 

call  dbinfo( "Define  Data=Surface  Centroids  List_Type=Siirface_List") 

Each  surface  in  the  list  has  three  real  values,  a  location  in  three  space  or  an  offset 
from  one  of  the  comers.  It  depends  on  how  the  data  is  used  by  the  application  pro¬ 
gram.  The  data  base  does  not  know  the  difference. 

call  dbinfoC'Define  Data=Surface_Centroids  Data_Dim=3  valuejype=REAL”) 

All  three  of  these  commands  could  be  combined  into  a  single  command.  (Pretend  all 
of  the  keywords  are  in  one  string  and  do  not  have  embedded  c;irriage  returns.) 

call  dbinfoC’Define  Data=Surface_Centroids  Data_Type=LlST 
List_Type=  Surface_Li.st  Data_Dim=3  value_type=REAL’’) 

Later  when  we  find  out  how  many  surfaces  there  are  on  this  object,  the  Surface_List 
can  be  defined.  If  there  are  1,001  surfaces,  then  the  data  definition  can  be  completed. 

call  dbinfoC'Define  List  Def-Surface  List  List_Length=100l") 
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3.2.4.  Example:  Definition  of  Grid  Dependent  Data 


The  Task 

Define  a  grid  dependent  data  item,  "lon  density",  which  is  element-centered. 

The  Solution 

As  in  the  example  before,  the  data  item  definition  is  independent  of  the  grid  and 
element  definitions.  The  data  is  defined  as  being  of  data  type,  "SPATIAL".  This 
implies  a  different  set  of  data  is  required  for  each  grid  in  the  problem.  SPATIAL  data 
is  defined  on  the  level  of  an  element.  Grids  are  thought  to  be  rectangular  lattices  with 
equal  spacing  in  each  direction.  An  element  is  located  at  each  lanice  point.  One  of 
the  predefined  data  items  is  one  called  "ELEMENT_CENTERED'’.  We  can  save  some 
work  defining  Ion_density  by  taking  advantage  of  an  existing  definition  using  the 
"SAME_AS  keyword". 

call  dbinfoC’DEFINE  Daia=lon_density  Same_As=ELEMENT_CENTERED") 

The  element  type,  ELEMENT_CENTER,  is  also  predefined. 

3.2.5.  Sample  Routine  1 

This  routine  demonsuates  how  to  use  statement  addressing  functions  to  use  the 
data  in  the  data  base  manager. 
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■:::v:ra::‘:-r  "SO  scar 


■V  ..-J  f 


‘ ,  I  z  )  -  L  v' :  o  -■ " 


w’  I  b  p ; 

( iz-  1  ) 


: ) 

» 


rigrids  =  idir^gd 
i)Cp  on  girio.3 

do  700  iorid  =  ] 


^  ^  ^  o 


cc 


:c  *reaa  in  pctenoia-s 
2000 


v;r  ite  (scorrnd,  2000  )  inr  i  0 

forn.ar.  {  'READ  DATA=p':.-._2''_ncde  gr :  d-^ '  ,  i  >-' ; 
call  dbaata(  sconn  i  ) 


'grab  stuff  out  of  output 
iofset  =  idboff(l) 


r.  (cdbdata.'r.i  wr.i;!.  will 


c 

2500 


get  inforo’.ation  about  tr.o-  uo cent ial s 
write  ( scorr.nd,  2  50  C  /  iorid 


format ( 


'  T  M.'"'*’’  ''  13  TT 


AATA-cot  20  nod€ 


r  io) 


call  dbinfo(  sccmnd  ) 
nxpct  =  idihir(l) 
nypot  =  idihir(2) 
ntpot  =  idihir ( 3 ) 

c  *locp  on  nodes  of  potenoials 
do  600  iz=l,nzpct 

600  continue 

c  *all  done  with  this  grid,  write  : 

write  (  scomnd,  7  000 )  igrid 


out  and  release  its  allocat 


iOO  format  (  'WRI' 

call  dbdata( 
;0  continue 


CLOSE  DATA=pot_20_nc de  grid=',i6) 
:omnd  )  "  ” 


'  Cj  f  lit"' 


enc 
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10.5  SAMPLE  ROUTINE  2 


This  routine  demonstrates  how  to  write  a  wrapper  routine  to  hide  the  data  base 
completely  from  the  calculation  routine.  foobar{). 


•r 


ne  local  variacles 
character’s^  sttr-.cd 


cc  ’t.tr.d  out  about  grids,  save  the  number  of  grids 
call  dbdata (" Inquire  3rid=ALL") 

cc  ’loop  or.  grids 

do  "^00  igrid  =  1,  r.grids 

cc  'read  ir  pctentials  for  this  grid 
write  (sccrtnd,  2000)  igrid 

20C0  format!  '?ZP.O  DATA=pot_20_r.c. de  grid=',i6) 

call  dbdatal  sco.mnd  ) 

c  ’grab  stuff  out  of  output  commo.r  (cdbdita.h)  which  will  be  useful  1 
iof set  =  idbcf  f ( 1 ) 

cc  ’get  informaticn  about  the  potentials  in  this  grid 
write (sccmrd, 250C)  igrid 

25CC  format!  'INQUIRE  DATA=pot_2C__node  grid“',i6) 

call  dbi.rfc!  sccmrd  ) 
nxpot  =  idih ir ( 1 ) 

rypot  =  idihir(2) 

rzpot  =  idihirO) 

c:  ’call  the  calcula’ion  submodule 

call  f oobar ( r data ( lof set ) ,  nxpot, nypot , ncpot ) 

cc  ’all  done  with  this  grid,  write  it  out  and  release  its  allocated  m.e 
write ( scomnd, "000 )  igrid 

7CC0  format!  'WRITE  CI-OSE  DATA--pot_20_node  grid=',i6) 

call  dbdata!  sccmrd  ) 

711  ccrtinue 

return 

end 

subroutine  focbar (pct20s, nr , ny, nz) 
real  pot2Cs(  nx,  ny,  nz  ) 

return 

end 
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Data  ite"  r.a.r.es  r.ay  h-.ve  multiple  flelus.  The  r.ar-e  f:r  a  ciete 
cf  data  will  be  the  data  type  name,  ftllcwed  by  the  rest  the 
fields  separated  by  one  space.  This  constructed  name  will  be 
truncated  to  80  characters.  If  truncation  cf  a  name  occurs,  an 
error  message  is  generated.  When  defining  a  data  tte.m.  wtth 
multiple  na.m.es,  just  use  the  first  nam.e  to  set  its  prccerttes. 

Then  when  making  dbdataO  calls  cr  inquires,  use  the  entire  nam.e. 

For  now,  spatial  ite.mc  are  3  dimensional.  When  the  keyword 
reader  gets  smarter,  keywords  using  input  lists  will  u.5e  (  ,■  to 
enclose  the  list.  For  example,  "FT ';£_LF;j3TH=  ( len_:-; ,  len_y ,  len_c  )  "  . 

Output  is  returned  via  a  include  file.  It  is  unwise  to  rely  cn 
the  order  cf  the  contents  of  the  include  files  since  tney  m.ay 
change.  The  names  cf  the  variables  in,  the  file  should  be 
constant . 

It  should  be  noted  that  the  current  implem.entation  cf  ti.m.e  names 
has  been  done  in  such  a  way  that  it  m.ay  be  a  useful  guide  for 
future  modif  icat  ions .  Timie  sta.mps  are  sim.ply  nam.es  which  can 
be  used  to  group  otherwise  identically  identified  data  item^s . 

If  another  cr  additional  distinguishing  r.am,es  were  desired, 
the  t  im.e  sta.m.p / namie  usage  could  be  generalised  t  o  all  ow  a 
nu;nber  of  different  and/or  parallel  sets  of  data  grou 
The  relationships  of  the  groupings  would  be  determ.ine 
a  set  cf  submodules  or  functional  re*l  at  ionships . 


BUGS:  All  err.crs  are  fatal  right  now. 

U  c  t  a  1 1  c  cmm.a  n  d  s  are  fully  i  mp  1  eme  n  t  e  d  . 
D e a  u  1 1  s  a r  e  n c  t  a  1  w a  y s  dec ume n  t e d  . 
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U.  -O 


:llowing  notes  are  sob;:ect  to  char.je  aj  aev< 
should  be  faurl/  ttable. 


ocr.t  i.nues, 


:he  kevucrds  have  been  mar<ed  to  indio:;*e  implementation  status. 
Mo  mark  means  it  is  i-mplemente-ii  and  may  even  wcr:-i. 

;  -  Partially  i n.p '.emented .  Restrictions  app.  y. 

■<-  -  Partially  irrn  1  emented .  Reooanited,  but  lur.ored. 


.\c 


irr.ple.~e.ntei.  Will  generate  fatal  er: 
;3  are  ir.  ;  ]'s,  where  appropriate. 

.u»  o  ^  t"  h'  ^ 


J  ,  ,  *  C.  wv  -tl 


1  -e.Tb 


-  *1 1  r.  ^  s 

^  V-  p:,  - 


described  at  t: 

no  *"  h p  data  r  'j 
on  s,3dawn  in 


.e  (; 


;....s  are  ava:..3D^e  in  tne 
.Dac/Motes/ruamcles . 


l.ne  ter  n-e  r.a:: 
.ented  inou 


O' _ -  . 


'.d)  -  i.nterface  r: 
rna  is  a  keyword  : 

'Ipen  ?refix=Dmsc" 

'STATIC  pref ix=Dat a_3ei  _  _  _ 

ru*:  iron,  the  routi  ne  is  returned  ir.  "c'ifcf  ile  .h" 


*i.<e: 

1  PREFlX^Data  set  2" 


ranu  xeyweros  are  ; 


EXIT  - 

Close  all  open  files. 

OPEN  - 

Assign  file  prefix  (s) 

to  run 

-  Close  file  prefix(sj 

. 

STATUS 

-  Query  current  file 

status 

♦ 

WHAT  - 

ST.ATUC  of  all  known  o 

refixes 

♦ 

HELP  - 

This  list. 

ifying  in 

f : 1 mat ion : 

PREFIX 

=File  Name  -  Define  a 

file  pr 

FILE  T 

Y?E‘‘keyword  -  Specify 

a  file 

SAVE  F 

ILE^logical  -  save  fil 

e  when 

* 

DIAGNC 

STIC=keyword  -  turn  ON 

, [OFF i 

ne  pretix.  [^asc  preri 
ype .  (DYNAPAC) 


X  u  s  e  a : 
ete (F a: 


;f) 


:a ( cccmnd) 

p  rp  ^ 

"rea 


.  s  a 


lerface  routine  for  date 
keyword  oriented  input 


base 

string 


a  data=pot_2  0_ 
from  the  rou^u 
3  are  broken  : 


q  >-  T  ri  =  0  " 


output 
key wc  r : 

1  instructions: 

OPEN  -  .Allocate  space  fc: 
CLOSE  -  Teallocate  soace 


:ia= 

retu: 


ned 


requests 

like 

)dbdata . h" 


jups  below; 


data  in  memo r y . 

;or  data  in  memory. 


CRUNCH  -  Squeeze  the  q.3ps  out  of  memory. 
READ  -  Transfer  data  from  disk  to  memcry. 

( '■  -  must  open  disk  prefi:-:  fir.ot  ) 
WRITE  -  Transfer  data  froTi  memory  to  disk. 

( -  .must  open  disk  prefix  first) 
COPY  -  o-ny  first  data  name  t- 


est  of  data  names  (in 


-ry)  . 
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cpiiiona'  ccn'jT'.ards  to  avemae  data  base  i~t  alt  actitr.s; 

lENd-TH  =  ir.t  -  Site  ir.  wards  r.eedei  rar  data  ir.  mer.trv. 


to  rr.cdify  data  type  def  rr.rr:  icr.  terr.porar  ily  .  :::r~ally 
this  keyword  is  only  used  ta  ir£N'  data  iter.s  wrtr. 
rr.ultiple  field  names.  In  •hir  rase,  the  data  type 
is  defined  with  by  the  first  field  only  while  Icraticns 
in  memory  are  associated  with  the  full  name. 

.HERE  =  int  -  Point  to  a  memory  address,  use  idblrcO  tr  pet 

offset  from  idata(CK  Used  to  tell  the  data  base  where 
something  is.  Note  that  m.urh  of  the  internal  error 
cherking  will  be  defeated  by  this  command. 


Grid  and  Data  Def initicn/Type  Oriented 

dbinf o  ( ccomr.d)  -  interface  routine  for  defining  and  learning 
information  about  grid  definitions  and  data  typin 
properties.  Sorry,  only  one  PROBLEM,  PRID,  EDEME 
or  INTERP  per  call.  Use  dbdata  ( "£T.'-.I'U3  DATA=stuf 
to  fi.nd  the  memory  offset  for  the  "stuff"  data  ir.  grid  2, 

a  dbinfo()  ccom.nd  is  a  keyworsd  cri-vr.tea  input  string  like. 

"DEFINE  Data=pot_8^node  elem.ent_type=Ncde_8 

dimiension=3caiar  va lu-'-real  set_value=C.  " 

output  from  the  routine  is  returned  in  "odbinfo.h" 

command  keywords  are  : 

DEFINE  -  Define  or  redefine  a  piece  (s)  of  inf or.maticn . 

INQUIRE  -  Request  information  about  a  GRID,  ELEMENT,  D.AT.d,  cr 
LI5T_DEF,  INTERP.  (see  context  keywords* 
DIAGNCSTIC=keywcrd  -  turn  ON,  iOF.-';  d_agncstics 
*  HELP  -  This  list. 

context  keywor'ir, : 

PROBLEM  -  Get  the  gene:  illy  userul  stuff  about  proble-  . 
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1 1,  V:  kQ 


CO  o 


DATA~n3.rr.c  ~  C~T^.3.r.cis  3.*^  '  '^.  is  ri^**  3  *'yr  'r 


INT:LRP=narr.e  -  net  implement^a  .  "  linear"  ,  "quaarae  _c" 
3Ai^’E_AS  =  na.T.e  -  T'ne  new  thini  is  tne  "same  as"  name. 


IS 

f  elements  in  a  r 
aa  describing  wh,' 


I ram.eter s  : 

DATA  TYTE’=k.evwcrd  -  tnis  da' 

CATA_LE1^GTH- inteaer  -  numbe; 

ELEMENT_Ty?  E  =  e  1  ement_narr.e  - 
l;st_TY?E=1 i3t_name  -  name  c 
D  AT  A  OX  M~  O "  r”  ^  '  3  C5  ■“  0  i.  t.  *Pi0  3  '"•  Mrr  ir  ^  ( V3 ' '  s  /  **  )  ^  ^  'j 

VALUE  TYPE—val  :e  t:yoe  ~  tyoe  cf  values  at  each  point 

V  a  u  e  t  y  ^  e  ^  R  e,  a_;  /  4.  ^  ^  e.  r*-, ,  l.  vyO  i.  u  A  —  /  ..  .  P*.  ^  /  n  ^  P  "  A.  ^ 

AL?HA_LE:JGTH  =  rchar3  -  length  c:  ATPKA  types  in  characters  [SC 
SET_VALUE=va l-.;e  -  def  ault  in i ••  i  a  1  i rat  icn  va  lue/.m-athod  . 

( ''  "  cn^. y  Kno'w  hca-  t  c  se*?  ^A..e _ 1  *!<■=...  t  '  censtants} 

A’lE_DATA=l':gi -al  -  save  data  .-.n  per.m .  [  TREE ’  cr  te.mp.  (FALSE) 
w;':_riLE  =  keyA- r d  -  [EA^SE[  -  put  this  data  in  the  main  f'^le 

-  TRUE  -  c  am-;  as  D':r;A?AC  cel  aw 

-  DYGAPAC  -  cut  data  in  own  s-.andard  subfile 

-  L0NG_NUM3ER  '':vimber_subkey s>  -  for  data 

types  with  many  subkeys  (  >'20C) . 

<nuntcer  G'..;b<ey3>  is  m.axirr.um  number  of 
k 6 y s  1 2*  cut  in  t h 0  f  i  1 0  ( d 0 f  0 u i t  is 
ICOC).  Addressable  only  with 
DATA_MAi'^E  <numter>,  where  <num.ber'»  is 
between  1  and  <nu.m.ber  subkey3>. 
FllE_SUFFIX  =  n.ame  -  name  to  use  for  OWN  FILE  suffi;-;  (should 

include  "."  as  first 
DRID_DEPEND=logical  -  TRUE  if  data  is 


T :me_depend  =  1 1  g 1 ca : 

TIME  SAVE  =  inteofer  - 


-  TRUE  i: 
hew  many 


Jata  IS 


i.T.e  steps 


character  ) 
grid  dependent  . 
time  decendent . 
data  ■ 


I 


Automat ic ally  OVER_W'RITEs 
versiin  of  data. 


Oldest  (smallest  t 


save  [1] 
ime  stam.t 


'r...  oarameters: 

UN:T_SIZE=integer  -  how  many 
UNIT_OFFSET  =  o  f:.3et_x,  c£fse':._' 
from  ele.mer.- 
d  1  m  e  n  5  _  c  n  a : . 
e.ach  "UN:T_.C 
rNTERP_TYPE  =  name  -  me*:hoc  of 

r_LEF  parameters : 

LIST_L.FNGTH  =  intege:  -  nir^ber 

iR?  parameters:  nene  for  now,  ;u. 

rr.aybe  later  will  id  a:-c’"^3+r. 


'/  a  ^  u  e  s  a : 


r :  3  0 1 


eacr.  unu- 
mesh  unit 
appr  Dp . 
offset  fot 


r  r  e  s  r.  c  u ..  u  i'  e  o  ri  • 
ftr  SPAT  (A!  “  ems  . 

;rDo:a*:icn  between  ccints 


entries  in  tne  ..i3t_nam,e 
the  na.me .  Will  recognize 
1 .  .  .  t  yce  Gt  u i  :  . 
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MJ 


GRID  parameters  {note  GPTD  :s  i  lis"  at  -'PA^TD',  elements); 

OUTSIDE  -  this  grid  en  :ltses  th<r  other  grids  (if  any 
INS  ICS  =  gr  id^n.ame  -  this  gria  is  ir.siie  srij_r.i;''e 
GRID_SIDE=reaI  -  gr_d  mesh  site  (meters),  redefines  t-l  ' 
MESH__RATIG=ir.teger  -  ratio  of  inner /cuter  mesn  units  ;>=G 
PRiM__RATIO=intecer  -  ratio  cf  i..ner /primary  (cr  cuter  ms; 

iTi©  3  Ti  U  P.  2- 1.  S 

ORIGIN=vaiue_x,  vaiue__y ,  value_z  -  In  the  INSIDE  grid  tne 

origin  is  in  cutergrid  units.  For  CISll 
nrids,  tne  origin  is  the  location  of  tne 
grid  1  in  this  grid. 

(  always  true;  CRIG  :i:  =  l ,  1 ,  1  in  cuterm.cst 
ORIGIM_?R IM=val’ue_:-: ,  va  1  ue_y ,  va lue_z  -  Drigi.n  Itcatiin  i.n 
pri.mary  cr  cuter  .most  grid  u.nits. 
EDGZ_LENGI.H=length_;-:,  length_y,  lencr  h_z  -  numiser  of  node  ; 

along  the  grid  edge,  includes  endpoint-, 
i  GRID_INTERP  =  na.m!e  -  cuter  grid  interpclatio.n  funottc.-. 

(1  -  =0  for  new) 


Th6 

following  items  are 

predefined  : 

(  < 

ntne>  means  no  defau 

It  value 

is  set 

for  that  att r i'e 

ute  ) 

Cat 

a  Types  [OwnFile=FAl 
BUFFER  Data  types 

SE] 

Name 

DatTyp 

VaiType 

Len  SaveDat  I 

nitVal 

Dim  GrdDep 

3UFFER_DATA 

BUFFER 

INTEGER 

C  fal.se  ' 

0' 

1  FALSE 

LIST  Data  types 

Nam.e 

DatTyp 

ListTyp 

VaiType  SaveDa 

t  InitV 

a.  -  D 1! 

LIST_DATA 

LIST 

<none> 

INTEGER  TRUE 

'  0' 

1 

F.AL 

SP.ATI.AL  Data  types  (ValType= 

REAL,  In 

itVal='  0.0',  S 

a veDat= 

TRUE , 

Data 

Gr idDep=TRUE, 

TimeDep=FALSE, 

NumTimSav=l j 

Data  Type  Name 

Element 

Type  Name 

SPATIAL  DATA 

<nor.e> 

NODE  8 

NODE  3 

ELEMENT 

NODE  20 

NODE  20 

ELEMENT 

NODE  32 

NODE  32 

ELEMENT 

ELEMENT_CENTERE 

D 

ELEMENT 

_CENTER 

£2.6 

.ment  Types 

Name 

Pts/Elm 

Interp 

Pt  Logs 

ELEMENT  CENTER 

1 

LINEAR 

(0.5, 

0.5, 

0.5) 

NODE  8  ELEMENT 

1 

LINEAR 

(G.O, 

0.0, 

0.0) 

NODE  20  ELEMENT 

4 

QUADRATIC 

(G.O, 

0.0, 

0.0) 

(0.5, 

0.0, 

0.0) 

(0.0, 

0.5, 

0.0) 

(0.0, 

0.0, 

0.5) 

NODE  32  ELEMENT 

4 

QUADRATIC 

(0.0, 

0.0, 

0.0) 

(0.5, 

G.O, 

0.0) 

(0.0, 

0.5, 

0.0) 

(0.0, 

0.0, 

0.5) 

List  Types  (Le.ngths  should  be  redefined) 


i.mDep 
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(/)  CO 


^rT'.0 

'JRFAC£_LIST 
'JR?ACE_NODES 

(the  outermost  grid  should  be  redefined) 

Name  MeshLen  Origin  EdgeLen  Grid_Inter 

'i'  1.  m  (1.,  1.,  1.)  (7,  9,  11)  NONE 


file  "DBhead" 


Length 
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3.3.  Modifyinn  Screens 

The  ceding  lor  the  screens  exists  as  ascii  tiles  in  SI)'t’NAPAC/inputs/*/*.sc  It  is 
possible  to  niodit’y  these.  There  is  currcnily  no  documentation  on  this  process 

3.4.  Potential  Interpolation  Schente 
3.4.1.  Continuous  Field  Interpolants 

We  originally  proposed  using  triquadratic  interpolants  in  DynaPAC  in  order  to 
obtain  more  nearly  continuous  electric  fields  than  were  provided  by  the  trilinear  inter¬ 
polants  used  in  previous  three  dimensional  codes.  Experience  has  shown  that  the  pro¬ 
posed  method  sutters  from  the  two  disadvantages  that  ( I )  the  fields  are  still  not  strictly 
continuous;  and  (2)  the  dilTerence  in  weight  between  the  corner  nodes  and  the  edge- 
center  nodes  leads  to  a  poorly  conditioned  matrix  Therefore,  we  have  abandoned  the 
triquadratic  formulation  in  favor  of  a  finite  element  interpolation  scheme  which  has 
strictly  continuous  electric  fields. 

The  basic  interpolation  functions  consist  of  functions  which  have  unit  value  and 
zero  slope  at  their  home  nodes,  and  zero  value  and  slope  at  opposite  nodes: 

F,,  =  (/-!)'/( I -2z-t-2z') 

F,  =  (z)7(  l-2z-i-2z‘) 

and  functions  which  have  zero  value  and  unit  slope  at  their  home  nodes  and  zero  value 
and  slope  and  opposite  nodes: 

C,,  =  z(l-/)F„ 

G|  =  /.{  /.- 1  )F, 

The  functions  F^  and  G,,  are  shown  m  figure  2.6.  It  can  be  verified  that  these  func¬ 
tions  can  exactly  reproduce  constant,  linear,  and  quadratic  potentials. 

We  generalize  to  three  dimensions  by  assigning  to  each  node  of  the  unit  cube 
four  interpolation  functions  corresponding  to  the  potential,  x-component  of  potential 
gradient,  y-component  of  potential  gradient,  and  z-component  of  potential  gradient  for 
that  node.  Thus,  the  contribution  of  these  four  quantities  at  point  fO.O.O]  to  the  poten¬ 
tial  at  [x,y,z|  is  governed  by  the  interpolation  functions  F,,(x)F,,(y)F,|(z). 
G„(x)Fj)(y)F„(z).  F,,(x)G„(y)F„(z).  and  F„{x)F,|(y)G, ,(/.).  (Contributions  of  other  cube 
comers  are  obtained  by  changing  "0"  subscripts  to  "I"  as  appropriate  )  We  elect  to 
omit  terms  of  the  form  GFG  or  GGG. 


Interpolation  Functions 


Dist.  from  Node 


Figure  2.6.  Interpolation  functions  Fq  and  Go- 
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3.4,2.  Interpolating  Fotentials  and  Fields  for  Special  Klenients 

Special  cleincnts  arc  those  clctiicnts  which  contain  siirtaccs.  Because  surface 
containing  cleincnts  arc  not  empty  cubes,  the  potential  interpolation  lunctions 
tleseribed  in  the  previous  section  cannot  lx:  straight Ibiwardly  applied.  !o  de^elop 
finite  element  matrices  tor  these  cells,  we  require  a  prescription  lor  calculating  the 
potential  and  electric  field  in  the  cell  interior.  Note  that  a  "special  element'  is 
bounded  by  three  types  of  surfaces;  (I)  square  surfaces  bounded  by  grid  edges  and 
shared  with  adjacent  (presumably  empty)  elements;  (2)  object  surfaces;  and  (3)  sur¬ 
faces  bounded  by  both  object  points  or  edges  and  grid  pt.'ints  or  edges.  On  type  ( I ) 
surfaces  we  must  use  the  potential  and  field  interpolation  described  in  the  previous  sec 
tion.  On  object  surfaces,  the  potential  and  electric  field  must  he  e.vpressed  in  terms  of 
the  object’s  surface  cell  potentials  and  normal  fields.  Type  1.3)  surfaces  must  smoothly 
blend  between  the  two.  We  describe  below  an  algorithm  to  express  the  potential  and 
field  at  any  point  in  the  cell  volume  in  terms  of  the  grid  point  potentials  and  fields, 
and  the  surface  cell  potentials  ami  normal  fields.  The  algorithm  has  the  propeny  that, 
when  applied  to  a  cell  with  six  type  ll)  surfaces,  constant,  linear,  and  quadratic  poten¬ 
tials  can  be  represented  exactly 

Let  R  be  a  point  in  a  volume  bounded  by  a  set  of  surfaces  j.S).  each  (d  which 
may  be  a  triangle  or  a  planar  convex  quadrilateral.  For  a  given  .S.  let  I*  be  the  point 
on  5  nearest  R.  let  ([>(1’)  and  K(P)  be  the  potential  and  electric  field  at  P.  and  n  fx;  the 
unit  normal  (from  the  surface  point  P  into  the  volume)  to  the  surface.  Let  At.V)  be  the 
area  of  the  surface.  Let 

d  =  R  -  P 
=  d  •  d 

K(.V)  =  A(.V)/  <r  for  d  •  n  >  0 

A  =  K(.V) 

Let  /  be  the  distance  from  P  in  direction  n  to  the  next  surface  intersection,  (In  prac¬ 
tice.  it  is  adequate  to  extend  /  to  the  intersection  with  the  surface  of  the  rectangular 
parallelepiped  bounding  the  volume  )  Then  the  contribution  of  .S'  to  the  potential  at  R 
is 

(j)^.(R)  =  (K(S)/N)  |(|)iP)  (I  -  d  •  n//)  (d  •  n)  (KtP)  •  n)] 
and  the  total  potential  is  found  by  summing  the  contributions  from  all  the  surfaces. 

The  electric  field  is  found  by  differentiating  thv  above.  The  contribution  of  ,S  to 
the  electric  field  is 

=  4)^.  (VWA  -  VK(.S)/K(.S')1 
+  (K(,S)//V)  I  V(t)(P)  +  1 1  2  d  •  It//)  n  ( KfP)  •  n)j 

where  V(|)(P)  is  taken  tangential  to  the  Mirlace.  so  that 
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V(|)(P)  =  0  tor  P  at  a  vciicx  of  S 

V<5)(P)  =  -e  (f  •  K(P))  tor  P  on  edge  with  direction  e 

V4)(P)  =;  -K(P)  +  n  (E(P)  •  n)  tor  P  in  surface  interior 


To  implement  this  we  interpolate  surface  fields  and  potentials  on  bounding  sur¬ 
faces  in  a  manner  which  maintains  the  continuous  lield  property.  Potentials  and  nor 
mal  field  components  are  defined  at  the  centroids  of  the  original  ohjeci  surfaces  and 
area-weighted  averaged  at  surface  corners.  Potentials,  normal  fields,  and  tangential 
field  components  at  all  surface  points  which  are  vertices  of  hounding  surfaces  arc 
expressed  as  linear  combinations  of  centroid  potentials  and  fields.  This  transformation 
is  given  by  the  matrix  TNdC'nt  in  the  special  element's  "Pound  Suils"  record. 

The  format  of  the  "Bound  Surfs"  record  is: 


Word 

Variable 

Contents 

1 

NC'Surf 

.Number  i  f  bounding  polygons 

2 

NCNoile 

Number  of  nodesla) 

.3 

N'C'nts 

.Number  of  centroids 

4 

1, Trans 

Length  of  constraint  relation  matrix 

next  8*NC.Surf 

JCSurftH.T) 

Bounding  polygon  node  list 

next  .'^^NCNchIc 

RCNodd.T*) 

Node  k)cations 

next  NCnts 

MatCntf*) 

Centroid  list 

next  LTrans 

I.TCntfM 

Consfrainf  relation  map 

next  LTrans 

TNdCntf*) 

('onsiraint  relation  matrix 

(a)  Nodes  l-K  are  cube  corner  nodes.  Remaining  nodes  are  surface  cell  centroids. 

surface  cell  corners  (constrained),  and  additional  surface  nodes  (constiamcd), 

3.5.  Surface  Potential  Boundary  ('onditions 

Physically,  it  is  possible  to  specify,  for  each  surface  element,  either  its  potential 
or  its  surface  normal  electric  field.  (For  conductors,  the  choice  is  between  potential 
and  total  charge.)  From  the  point  of  view  of  the  inteipolation  functions  (and  therefore 
the  potential  equations),  however,  the  potential  and  surface  field  ap|X’ar  as  independent 
vanables.  DynaPAC  legislates  the  physical  constraint  that  one  or  the  other  must  he 
specified. 

Fxpenenee  h;is  shown  that  the  surface  electric  fields  are  weakly  coupled  to  the 
potential  equations,  so  that  solving  for  the  surface  electric  field  variable  docs  not  give 
reliable  resulfs.  Rather,  wc  always  solve  Poisson's  equation  with  all  surftice  electric 
fields  fixed.  To  determine  surface  field  values  for  fixed  potential  cells,  we  use  the  fact 
that  the  surface  cell  charge  is  conjugate  to  its  potential,  i.e..  when  the  full  conjugate 
gradient  matrix  multiplication  is  performed,  the  surface  cell  charge  (corresponding  lo 
the  existing  potentials  and  fields)  appears  in  the  surface  cell-potential  slot  of  the  "AF" 


vector.  This  is  an  iniprovcincni  over  "eapaeitance  niairix  '  inelhods  in  tliai  i(  lakes  lull 
account  ot  the  nonlinear  space  charge  in  the  external  space 

The  matrix  operations  refened  to  alxive  are  earned  out  hy  Subrouline  1  (  optd 
The  surface  ceil  charge  values  are  then  passed  to  Subrouline  IT'pdal  to  dclennme  neu 
surface  electric  ticid  values.  The  new  value  is  lirsl  determined  b\  dividing  the  charge 
by  the  cell  area.  Then.  FTlpdat  mixes  of  the  previous  (old)  values  with  ihe  ncwlv 
calculated  values.  F  inally,  FiUpdal  bounds  the  allowable  electric  fields  bresed  on  the 
cell  dimensions,  the  minimum  mesh  si/c  in  the  problem,  and  the  range  of  potentials  in 
the  problem.  A  diagnostic  is  printed  for  those  cells  whose  calculated  electric  fields  fall 
outside  the  allowable  range,  and  the  unpei  or  lower  bound  (vis  appropriate)  is  iiscii  for 
the  new  electric  (leld  value 

Most  of  (his  machinery  appears  to  work  l)vnal’kl7ll'S  may  be  used  to  set 
potential  values  on  conductors,  and  potential  or  electric  licld  values  on  insulatini!  sur 
laces.  DynaF’Rh/llfS  may  also  be  user!  to  tieline  a  lloating.  biased  set  of  cimductors 
nutubered  consecutively  starting  with  ”1".  with  fixed  potentials  on  higher  numbered 
conductors.  Nonzero  charge  values  lor  these  eonduclois  may  be  generated  using 
fracker.  or  entered  manually  using  a  stand  ahme  module 

3.6.  Space  Charge  I  reatments 

The  various  space  chiu-ge  treatments  available  in  the  DynaPAC  potential  stdver 
were  enumerated  in  .Section  2.4.2  Here  we  give  some  addition  details  for  each  of  the 
options.  These  Ireatments  are  impicmenterl  in  subroutines  RhoOlP  and  PSSern 
RhoOlT  returns  values  to  be  used  for  the  charge  density  .mil  ils  derivali'  e  with  respect 
to  potential  I’SScrn  apportions  the  charge  (density  relumcil  by  RhoOlT  limes  ihc  cic 
ment  volume)  to  the  nodes  ol  iliv  element  as  if  the  charge  was  concentrated  ai  a  point 
in  the  middle  ot  the  element  fhe  derivative  is  treated  as  clemeni-centered. 

3.6.1.  I.aplacian  Space  Charge 

The  Laplace  space  charge  option  solves  Laplace's  ccju.ilion. 

-V’(j)  =  0 

i.e..  charge  exists  only  on  object  surfaces  and  external  boundaries,  as  detennined  by 
the  boundary  conditions.  "Space  charge"  iterations  may  still  be  required,  however,  due 
to  the  treatment  of  surface  electric  lields 


3.6.2.  Linear  Space  Charge 

The  Linear  space  charge  option  solves  the  Helmholtz  or  Debye  Huckei  equation 

The  value  of  X  is  the  maximum  of  the  plasma  Debye  length  (divided  by  the  square 
root  of  density  depletion  due  to  geometric  shadowing),  and  the  local  mesh  spacing  di¬ 
vided  by  the  DEBLIM  parameter  (which  defaults  to  2). 

3.6J.  Nonlinear  Space  Charge 

Nonlinear  space  charge  is  calculated  using  the  NASCAP/LEO  formula. 

-V>  =  p/e„ 

p/^'o  =  ( I  +  I  4>/0  I  C((|),E))  /  ( I  -1-  (4ji)  '  1  (J)/0  i  ''“) 

C{(|).E)  =  I  e/0  I  ((K  ./rr  -  1 1 

(R^j^/r)'  =  2. 2d  I  Fa'^/0  |  '  I  0/0  ! 

p  =  space  charge  |coui-m  '| 
e,,  =  8.H54xl()  |farad-m  '| 

Afj  =  plasma  Debye  length  |m| 

0  =  plasma  temperature  jeVl 
($)  =  local  space  potential  |  volts  ( 

E  =  local  space  electric  field  (volls-m  'l 

where  the  last  equation  (analytic  focusing)  comes  from  lilting  a  fmiie-lemperature 
spherical  (Langmuir-Blodgetl)  sheath.  If  analytic  convergence  is  lurned  off  (by  the 
CONV  input  parameter)  C(0.E)  is  set  to  zero.  The  value  of  is  modified  in  accor 
dance  with  the  mesh  spacing  as  described  for  linear  screening  (above). 

3.6.4.  Frozen  Ion 

The  "Frozen  Ion"  treatment  is  intended  for  short  timescale  (typically  a  few  mi¬ 
croseconds  or  less)  problems  for  which  it  is  a  good  approximation  to  assume  that  elec¬ 
trons  are  in  barometric  equilibrium,  but  ions  have  nol  moved  Ihc  space  charge  func¬ 
tion  depends  on  the  mesh  ciependeni  potential,  0|,  which  salislies 

I  -  exp(0|/9)  -  -(/v/D)"  ((])|/0) 

where  D  is  the  local  mesh  spacing  divided  by  the  DEBLIM  parameter  ((j)j  is  zero  for 
D</v.)  The  space  charge,  p/r^j,  is  then  given  by 

(0j/D" H  I -exp((t)/(!),)) 

-0/D"  ()>(})>(t), 

-0|/D'  1-  (0/X‘')(exp(0|/0)  -exp(0/0)) 
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3.6.5.  Tracked  loii-s 


This  algorithm  is  used  for  timesealcs  on  which  it  is  practical  to  treat  ion  motion, 
but  electrons  are  considered  in  banmietric  equilibrium  The  ion  density  has  been 
stored  in  the  RHO  Ion  array  by  the  Tracker  module,  and  the  electron  charge  densiis 
must  be  added.  .Subroutine  RhoOtV  returns  the  electron  charge  density,  as 

-((t)+e(L/X)T)/L~  0>() 

-(&}.")  exp((J)X‘/0L‘')  0^1) 

L  =  max(X.I)) 


3.6.6.  Explicit  PIC 

For  this  option,  it  is  assumed  that  the  Tracker  module  has  stored  both  the  electron 
and  ion  charge  densities  in  the  RHO  Ion  array. 

3.6.7.  Sheath  Wake  Densities 

This  option  is  intended  for  problems  in  which  the  "01 "  density  array  has  been 
calculated  using  shadowing  by  a  positive  potential  contour  (rather  than  by  the  object). 
The  module  which  perfonns  this  calculation  is  not  completely  developed.  The  section 
of  RhoOfP  which  treats  this  case  is  intended  to  blend  between  regions  with  accelerated 
electrons,  quasi-neutral  regions,  anil  electron  rich  wake  regions.  We  do  not  attempt  to 
explain  the  logic  here,  but  refer  to  interested  reader  to  the  RhoOlP  subroutine.  (.Sha¬ 
dowing  is  not  appropriate  to  large  negative  sheaths.  A  variant  of  the  "Tracked  ion" 
method  should  lx:  developed  for  such  cases.) 

3.6.8.  Special  Space  Charge  Formulation 

The  option  is  available  to  call  a  user-written  routine.  RSpecl.  for  a  developmental 
space  charge  option.  One  array  which  might  be  used  in  such  an  option  is  RHO  Elec, 
the  element -centered  electron  charge  density,  calculated  by  the  Tracker  module  b\ 
tracking  full  electron  orbits  in  existing  potentials. 

3.7.  Selected  Subroutine  Descriptions 

There  are  no  subroutine  descriptions  supplied  at  this  time  Just  to  provide  some¬ 
thing.  here  is  a  ileseription  of  the  "element  table".  LTBL.  which  is  written  into  the  Da 
taBase  for  each  grid.  For  each  volume  element  the  I  .TBI,  entry  is 


<0 

Element 

is  tilled  or  belongs  to  child  grid 

0 

Element 

is  empty. 

>0 

Special  e 

lenient  index  number 

4.  DynaPAC  Example  -  Part  I 

This  example  originally  appeared  in  the  Interim  Report.  SSS-I)PR-90- 1 1973.  30  Sep¬ 
tember  1990. 
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4.1.  Introducliun  to  Using  DynaPAC 

This  chapter  takes  you  through  a  sample  calculation  using  many  of  the  DynaPAC 
modules.  In  this  example,  you  will  use  DynaPAC  to  calculate  the  potenual  about  a 
charged  cube.  The  calculation  prtK'eeds  as  follows; 

1.  Use  PATRAN  to  define  a  cube  with  a  side  of  20  centimeters.  You  will  also  create 
an  auxiliaiy  material  file  that  defines  the  surface  material  to  be  aluminum. 

2.  Use  GridTool  to  define  a  grid  about  the  cube  with  one  level  of  subdivision. 

3.  Use  PatDyn  to  initialize  the  DynaPAC  data  base  and  place  in  it  the  needed 
geometrical  information. 

4.  Use  the  “main”  module  to  establish  an  initial  potential  (IPS)  of  -  KXX)  volts  on  the 
cube  and  to  create  an  input  file  for  the  potential  solver.  Among  the  input 
specifications  are  the  plasma  parameters  for  which  you  will  u.se  a  density  of  1  ¥ 
101 1  m-3  and  a  temperature  of  0. 1  eV. 

5.  Run  the  potential  solver. 

6.  Use  the  data  scanner  to  print  the  potentials  and  electric  fields  on  the  cube  surface 
and  in  the  space  near  the  cube  and  to  plot  the  potentials  in  space. 

4.1.1  Using  PATRAN  to  Define  the  Cube. 

Figure  4.1.1  .shows  the  PATRAN  session  file  to  define  the  cube. 
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Lines  1-2 
Lines  3-6 
Lines  7-8 
Line  9 
Lines  10-16 
Line  17 
Lines  18-19 
Lines  20-21 
Lines  22-28 
Line  29 


Sign  on  to  PATRAN,  requesting  a  new  datafile. 

Create  a  cube  of  the  required  dimensions. 

Define  “PATCH ’s”  on  the  faces  of  the  cube. 

Create  an  element  of  material  1  on  each  PATCH. 

Merge  the  equivalent  nodes  and  compact  the  ntxie  numbers. 
Check  for  free  boundaries. 

Display  normals. 

Correct  normals  to  be  consistently  outward. 

Write  the  neutral  file. 

Exit  PATRAN. 


1 

GO 

2 

1 

3 

PRIM.CUBE.BRICK 

4 

1 

5 

0.2.0. 2, 0.2 

6 

5 

7 

PRIM.CUBE.EV 

8 

S 

9 

MESH, P1T6, QUAD/4/1, N.1/1/1./1,-1 

to 

EQUIVALENCE 

1 1 

N 

12 

2 

13 

1 

14 

Y 

15 

6 

16 

2 

17 

VBOUND 

18 

VNORM 

19 

1 

2C 

3 

21 

1 

22 

INTERFACE 

23 

1 

24 

1 

25 

1 

26 

CUBE  FOR  SAMPLE  PROBLEM 

27 

N 

28 

N 

29 

STOP 

Figure  4  1,1  PATRAN  session  file  to  define  the  cube. 
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Finally,  rename  the  neutral  file  with  the  command 


mv  patran.out.l  cubc.ncu 

and  create  the  material  file  “cube. mat”  containing  the  single  line 

1  Aluminum 


4.2  Using  (iHd  l  oul 

1.  Issue  the  "GridTool”  command  to  run  the  GridTool  module  of  DynaPAC. 

2.  Choose  the  "Files"  option  in  order  to  specify  the  file  prefix.  Change  it  from  “fort” 
to  “cube.” 

3  Accept  the  default  of  a  PATRAN  neutral  file  with  units  of  meters. 

4.  Choose  the  “Edit”  option  to  change  the  grid  parameters.  Stan  with  the  only 
existing  grid,  grid  #1,  and  choose  the  "Modify"  screen.  The  console  now  appears 
as  shown  in  figure  1 1.2. 

5.  Alter  the  grid  dimensions  to  1 3  x  1 3  x  13  and  the  mesh  size  to  0.6  meters,  and 
accept  those  changes. 

6.  Next,  choose  the  “Add"  option  from  the  pull-down  menu  in  order  to  place  a 
subdivided  grid  around  the  cube.  The  screen  now  appears  as  in  figure  1 1.3. 
Change  the  subdivision  ratio  to  3,  and  the  grid  limits  to  run  from  5  to  9  in  each 
direction.  Accept  the  grid. 

7.  Choose  the  “Plot”  option  from  the  pull-down  menu.  Using  the  default 
parameters,  make  and  show  the  plot,  which  appears  as  in  figure  1 1.4. 

8.  Exit  the  “Edit”  menu,  select  the  "Write"  option  to  write  the  grid  definition  file, 
“cube.grd,”  and  exit  GridTool. 
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editor 
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Figure  4.2  GridTool  screen  for  the  Edit/Modify  option 


Show  Write  Reset  He 
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Figure  4.3  GndTool  screen  for  the  Edit/ Add  option. 


GridTool  plot  for  the  sample  problem. 


4.3  Running  PatDyn 


PatDyn  can  now  be  run  with  the  command 

PatDyn_SUN4  <  patdyn.in  >  patdyn.out  & 

The  first  line  of  the  input  file,  "patdyn.in",  uses  the  PREFIX  keyword  to  set  the  file 
prefix.  Normally,  no  other  lines  are  needed,  so,  in  this  case,  "patdyn.in"  appears  as 

PREFIX  cube 
END 

(The  input  file,  "patdyn.in",  can  also  be  built  by  the  DynaPre  module.)  The  output  file, 
"patdyn.out",  repeats  the  mesh  specification  information  and  contains  information  on  any 
possible  errors  that  may  have  been  detected  in  processing  the  object.  It  concludes  by 
giving  the  properties  of  the  materials,  although  these  are  not  currently  used  by  DynaPAC. 

4.4  Using  DynaPre  to  Prepare  to  Solve  For  Potentials 

Enter  the  DynaPre  module  by  issuing  the  command  “DynaPre”.  The  first  ta.sk  is  to 
specify  a  potential  of  -1()()0  volts  on  each  face  of  the  cube.  Choose  the  “IPS”  option  from 
the  top  menu,  and  choose  the  “Edit”  submenu. 

The  next  screen  offers  a  range  of  conductors  and  materials  to  which  initial  conditions  are 
to  be  applied.  By  default,  this  is  the  entire  range  available.  Accept  this  choice.  A  pull¬ 
down  menu  now  appears  offering  a  choice  of  the  type  of  boundary  condition  to  be 
applied.  Choose  “Fixed”.  The  screen  now  appears  as  shown  in  figure  4.5.  Change  the 
specified  potential  to  -KXX).  You  are  instantly  returned  to  the  main  IPS  menu.  Write  out 
the  initial  conditions,  and  exit  IPS. 
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Figure  4.5  DynaPre  screen  for  the  IPS/Fixed  option. 


Now  choose  the  “I^oient”  option  from  the  top  menu  and  the  “Hdit  Script”  option  from  the 
pull-down  menu.  You  will  be  presented  with  the  default  potential  solver  input  options,  as 
shown  in  figure  1 1.6.  Change  the  "Comment”  to  "Cube  Sample  Problem”  and  toggle  the 
"Problem  type"  to  "NON_LlNEAR”.  Note  that  the  density  and  temperature  are  already 
properly  set.  Exit  this  screen  and  choo.se  the  “Make  Script”  option  from  the  pull-down 
menu.  This  will  write  a  fully  commented  input  file  in  ”ps.Jn”.  Exit  the  pull-down  menu 
and  the  “DynaPre”  module,  and  inspect  this  file. 

4.5.  Running  the  Potential  Solver 

Run  the  potential  solver  using  the  command 

PotenLSUN4  <  ps_in  >  ps_out  & 

The  output  file  contains  a  reiteration  of  the  input,  problem,  and  grid  parameters,  and 
details  the  convergence  of  the  potential  solver. 

Before  proceeding,  note  the  “Fort”  and  "PS”  prefix  files.  "Fort.DP”  contains  object 
surface  information  and  various  miscellaneous  information  about  the  calculation. 
“Fort.BS”  contains  information  about  the  bounding  surfaces  of  the  special  elements  and 
is  needed  to  calculate  electric  fields  in  special  elements  or  to  plot  potentials  in  space. 
“Fort.ME02”  contains  the  potential  solver  matrices  for  special  elements  in  grid  2.  There 
will  be  one  such  file  for  each  grid  containing  special  elements.  "PS. DP”  contains 
information  generated  by  the  potential  .solver.  The  ‘  .til”  files  contain  the  data  base 
specifications  that  characterize  the  data  entities  in  the  other  files.  If  you  wish  to  save 
results  from  calculations  done  with  several  different  parameters,  you  need  save  only  the 
‘  .DP”  files;  the  remaining  files  are  not  changed. 
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Potential  Solver  Parameters 
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Figure  4.6  DynaPre  screen  for  the  POTENT/Edit  Script  option 


4.6  Using  the  Data  Scanner 


1.  Issue  the  command  “Scanner”  to  run  the  data  scanner.  Choose  “Print”  from  the 
top  menu,  and  “Data  Info”  from  the  pull-down  menu.  Change  the  "Working 
Grid”  to  2,  and  e.xit  the  “Data  Info”  screen. 

2.  Choose  “Edit  Prim”  from  the  pull-down  menu.  You  will  now  see  the  screen 
shown  in  figure  1 1.7.  You  have  the  option  of  printing  to  the  screen  or  to  a  named 
file  (which  can  later  be  edited).  Toggle  the  “Output  Mode”  to  “SCREEN".  Also, 
change  the  “Right”  and  “Bottom”  limits  to  9,  and  the  "Left”  and  “Top”  limits  to  3. 
Exit  the  screen. 

3.  Choose  “Make  Print”  from  the  pull-down  menu.  The  potentials  and  potential 
gradients  at  the  node  points  of  the  z=7  (middle)  gnd  plane  are  printed.  The 
potential  (volts)  for  each  point  is  at  the  intersection  of  the  appropriate  row  and 
column.  The  potential  gradient  (volt.s/meter)  in  the  y  direction  is  printed  above 
the  potential,  the  x-direction  gradient  is  printed  to  the  nght.  The  entries  at  (7.7 1 
are  all  zero,  as  this  point  is  inside  the  cube.  Note  that  the  potentials  and  fields  are 
consistent  and  that  there  is  a  slight  asymmetry  due  to  asymmetry  in  the  element 
bounding  surfaces. 

4.  Return  to  the  pull-down  menu  and  again  choose  “Data  Info”.  Specify 
“POT_Surf’  as  the  “Data  Name”.  Review  the  “Edit  Print”  screen  (it  should  be 
correct)  and  select  “Make  Print”.  The  potentials  (-1()(K)  volts)  and  normal  electric 
fields  (-5000  volts/meter)  will  appear  on  the  screen.  Again,  the  field  is  slightly 
asymmetric  due  to  asymmetry  in  the  bounding  surfaces.  Exit  the  “Print”  module, 

5.  Select  “Plot”  from  the  top  menu.  Review  the  “Data  Info”  screen,  and  choose 
“Edit  Plot”.  Figure  1 1.8  shows  the  “Edit  Plot”  screen. 

6.  Choose  “GRID  LIMITS”  and  change  the  limits  to  plot  from  3  to  9  in  each 
direction  (this  time  in  units  of  the  outermost  grid).  Exit  the  screen. 
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Figure  4.7  Scanner  screen  for  the  Pnnt/Edit  Print  option. 
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Figure  4.8  Scanner  screen  for  the  Plot/Edit  Plot  option 


7.  Choose  “CONTOUR  LEVELS”.  Change  the  number  of  contours  to  15.  The 
contour  levels  are  recalculated  instantly,  and  the  screen  is  updated  to  show  1 5 
levels.  Change  contour  levels  12-15  to  -30,  - 10,  -1,  and  0  volts,  respectively.  Exit 
the  screen. 

8.  Choose  “CONTOUR  MARKS”.  Change  “Number  of  Marks"  to  2.  The  .screen  is 
redrawn  to  show  2  marked  contours.  Change  the  second  marked  contour  value  to 
-30  volts  and  the  mark  to  “3”.  Exit  the  screen. 

9.  Choose  “LABELS  AND  INCLUDES”.  Change  the  “Title”  to  “Cube  Example”. 
Exit  the  screen. 

10.  Select  “OTHER  OPTIONS”.  If  you  are  on  a  color  terminal,  toggle  that  option  to 
“YES”. 

1 1 .  Return  to  the  pull-down  menu,  and  select  “Make  Plot”  and  “Show  Plot”.  A  gray¬ 
scale  version  of  the  contour  plot  appears  in  figure  1 1.9.  A  black-and-white 
version  of  the  contour  plot  appears  in  figure  11.10.  Hit  <RET>  to  go  back  to  the 
menu.  Exit  “Plot”  menu,  and  then  exit  Scanner. 

Note  the  presence  of  the  files  “fon.2”  and  “Scanner.opt”.  The  plot  data  is  contained  in 

"fort. 2”  and  can  be  displayed  or  printed  using  any  of  the  DynaPAC  graphical  interfaces. 

“Scanner.opt”  stores  the  options  used  when  ninning  the  data  scanner,  which  will  be 

restored  in  the  next  Scanner  run  in  this  directory. 

4.7  Summary 

This  example  has  taken  you  from  defining  an  object  through  calculating  and  displaying 

the  potentials  in  the  surrounding  space. 
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SIXV-A 
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Figure  4.9  Contour  plot  of  potentials  for  sample  problem,  gray-scale 
version. 
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5.  DynaPAC  Example  -  Part  11 

This  example  originally  appeared  in  the  Quarterly  Repon.  SSS-[)PR-9I- 12658.  July 
1991. 
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5.  A  Sample  Calculation  Using  D^naPAC 
5,1  General 

The  goal  of  this  calculation  is  to;  ( 1 )  calculate  the  biased  floating  potential  about  an 
object  consisting  of  two  disjoint,  unequal  size  cubes;  (2)  examine  some  particle 
trajectories  from  the  sheath  (ions  as  well  as  electrons). 

The  calculation  proceeds  as  follows; 

1 .  Use  PATRAN  to  define  two  cubes  having  sides  of  20  cm  and  30  cm,  respectively, 
and  being  20  cm  apan.  The  material  file  will  define  the  surface  materials  to  be 
aluminum  and  silver,  respectively. 

2.  Use  GridTool  to  define  two  nested  grids  about  the  object. 

3.  Use  PatDyn  to  initialize  the  DynaPAC  database  and  place  in  it  the  needed 
geometrical  infonnation. 

4.  Use  the  “DynaPre”  module  to  establish  a  biased  potential  set  where  the  small  cube’s 
potential  is  initialized  to  -10  volts  and  the  large  cube’s  potential  is  biased  -i-lO  volts 
relative  to  the  small  cube’s.  Also  create  an  input  file  for  the  potential  solver. 

5.  Run  the  potential  solver. 

6.  Use  the  data  scanner  to  print  the  potentials  and  electric  fields  on  the  object  surfaces 
and  in  the  space  near  the  cubes  and  to  plot  the  potentials  in  space. 

7.  Use  “DynaPre”  again  to  create  an  input  file  for  the  particle  generator.  Then  run 
p;uiicle  generator  (PartGcn)  to  generate  sheath  particles. 


84 


8.  Use  “DyiiaPre”  again  to  create  two  input  files  for  the  panicle  tracker  (Tracker).  One 
input  file  tells  Tracker  to  plot  the  particles.  The  other  input  file  tells  Tracker  to  track 
those  particles  and  generate  a  trajectory  plot.  Then  run  I  racker  twice  with 
appropriate  input  file  as  described  above. 

5.2  Parameter 

The  parameters  used  in  the  calculations  will  be: 


Pla.sma  density 

1  X  10'» 

m'^ 

Plasma  temperature 

0.1 

eV 

Plasma  species 

O+.e- 

Magnetic  field 

4  X  10-5 

tesla 

5.3  Using  PATRAN  to  define  the  cube. 

Refer  to  sample  1  for  how  to  define  an  object  using  PATRAN.  Define  the  object  (two 
cubes).  Finally,  at  the  UNIX  prompt,  rename  the  neutral  file  with  the  command 

%mv  patran.out.l  ZCubes.neu 

and  create  the  material  file  “ICubes.mat”  containing  two  lines 

1  Aluminum 

2  Silver 

5.4  Using  GridTool 

1 .  Issue  the  “GridTool”  command  to  run  the  GridTool  module  of  DynaPAC. 
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2.  Choose  the  “Files”  option  in  order  to  specify  the  file  prefix.  Change  it  from  “fort”  to 
“2Cubes.” 

3.  Accept  the  default  of  a  PATRAN  neutral  file  with  units  of  meters. 

4.  Choose  the  “Edit”  option  to  change  the  grid  parameters.  Start  with  the  only  existing 
grid,  grid  1,  and  choose  the  “Modify”  screen.  The  console  now  appears  as  shown  in 
figure  5.1. 

5.  Alter  the  grid  dimensions  to  13  x  13  x  13  and  the  mesh  size  to  0.16  meters  and 
accept  those  changes. 

6.  Next,  choose  the  “Add”  option  from  the  pull-down  menu  in  order  to  place  a 
subdivided  grid  around  the  two  cubes.  The  screen  now  appears  as  in  figure  5.2. 
Change  the  grid  limits  to  run  from  4  to  10  in  the  X  direction  and  from  5  to  9  in  the 
Y  and  Z  directions.  Accept  the  grid. 

7.  Choose  the  “Plot”  option  from  the  pull-down  menu.  Using  the  default  parameters, 
toggle  “PLOTTER”  to  your  favorite  plot-  reader.  Make  and  show  the  plot,  which 
appears  as  in  figure  5.3. 

8.  Exit  the  “Edit”  menu,  select  the  “Write”  option  to  write  the  grid  definition  file, 
“2Cubes.grd,”  and  exit  GridTool. 

5.5  Running  PatDyn 

PatDyn  can  now  be  run  with  the  command 

%PatDyn_SUN4  <  patdyn.in  >  patdyn.out  & 

The  first  line  of  the  input  file,  "patdyn.in,”  uses  the  PREFIX  keyword  to  set  the  file  prefix. 

Normally,  no  other  lines  are  needed,  so  in  this  case,  “patdyn.in”  appears  as 


PREFIX  2Cubes 
END 
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(The  input  file,  “patdyn.in,”  can  also  be  built  by  the  DynaPre  module.)  The  output  file, 
■‘patdyn.out,”  repeats  the  mesh  specification  information  and  contains  information  on  any 
possible  errors  that  may  have  been  detected  in  processing  the  object.  It  concludes  by 
giving  the  properties  of  the  materials. 

5.6  Using  DynaPre  to  prepare  to  solve  for  potentials 

Enter  the  DynaPre  module  by  issuing  the  command  “DynaPre.”  Enter  ■‘2Cubes”  at  the 
“File  Prefi.K;”  prompt.  The  task  is  to  specify  a  biased  floating  conductor  set  where 
conductor  1  (small  cube)  gets  initialized  to  -10  volts  and  conductor  2  (large  cube)  is 
biased  -rlO  volts  relative  to  conductor  1. 

Choose  the  “IPS”  option  from  the  top  menu;  then  choose  the  “Conductors”  submenu,  and 
choose  “Set”  option. 

The  next  screen  shows  a  list  of  all  conductors  defaulted  to  “Fixed  Potential.”  Toggle 
conductor  I  type  to  “Bia.sed  Potential”  and  set  its  “Value”  to  - 10  volts  (leaving  “Biased 
From”  to  be  “0”  since  this  is  the  lowest  conductor  number  in  the  set).  Toggle  conductor  2 
type  to  “Biased  Potential,”  set  its  “Value”  to  +10  volts,  and  set  "Biased  From”  to  “1." 

The  screen  now  appears  as  in  figure  5.4. 

Now,  go  back  to  the  main  IPS  menu  (by  hitting  <ESC>  twice).  Select  “Write”  to  write 
out  the  initial  conditions,  and  exit  IPS.  An  ascii  file  named  “ips.out”  is  also  written  out 
for  diagnostic  purposes. 

Now  choose  the  “Potent”  option  from  the  top  menu,  and  the  “Edit  Script”  option  from  the 
pull-down  menu.  You  will  be  presented  with  the  default  potential  solver  input  options,  as 
shown  in  figure  5.5.  Change  the  “Comment”  to  “Biasing  Cubes  Sample  Problem,”  and 
toggle  the  “Problem  type”  to  “NON  LINEAR.”  Note  that  the  density  and  temperature  are 
already  properly  set.  Exit  this  screen  and  choose  the  “Make  Script”  optio.:  from  thv.  pull¬ 
down  menu.  This  will  write  a  fully  commented  input  file  in  "ps_in.”  Exit  the  pull-down 
menu  and  the  “DynaPre”  module,  and  in.spect  this  file. 
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5.7  Running  the  potential  solver 


Run  the  potential  solver  using  the  command 
%  Potcnt_SUN4  <  ps„in  >  ps_out  & 

The  output  file  contains  a  reiteration  of  the  input,  problem,  and  grid  parameters,  and 
details  the  convergence  of  the  potential  solver. 

Before  proceeding,  note  the  2Cubes.*  files  created  so  far. 


2Cubes.DP 

Main  data  file  containing  most  of  information  about  the 

calculation. 

2Cubes.BS 

Contains  information  about  the  bounding  surfaces  of  the 
special  elements  and  is  needed  to  calculate  electric  fields 
in  special  elements  or  to  plot  potentials  in  space. 

2Cubes.ME02 

Contains  the  potential  solver  matrices  for  special  elements 
in  grid  no.  2.  There  will  be  one  such  file  for  each  grid 
containing  special  elements. 

2Cubes.HI 

Contains  the  database  specifications  that  characterize  the 

data  entities  in  the  other  files. 

If  you  wish  to  save  results  from  calculations  done  with  several  different  parameters,  you 
need  save  only  the  ".DP”  file— the  remaining  files  are  not  changed. 

5.8  Using  the  Data  Scanner 

1 .  Issue  the  command  "Scanner”  to  run  the  data  scanner.  Again,  enter  "2Cubes”  at  the 
“File  Prefix”  prompt. 

2.  Select  “Print”  from  the  top  menu  and  choo.se  “Data  Info.”  Change  "Data  Name”  to 
“POT_Surf  ’  and  exit  this  form. 
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3.  Select  “Edit  Print”  and  toggle  "Output  Mode”  from  “FILE”  to  “SCREEN  ”  Exit  this 
form  and  seleci  ‘Make  Print,”  The  screen  will  display  all  surface  potentials  and 
fields  as  in  figure  5.6.  Note  that  the  total  net  charge  on  the  two  cubes  is  almost  zero. 
Exit  this  screen  and  go  back  to  main  Scanner  menu. 

4.  Select  “Plot”  from  the  top  menu  and  choose  “Edit  Plot.”  Figure  5.7  shows  the  ‘Edit 
Plot”  screen. 

5.  Toggle  “Window  Manager”  to  the  most  appropriate  one.  Toggle  “Plotter”  to  your 
favorite  “fort.2  plot-reader.” 

6.  Choose  “GRID  LIMITS.”  Change  limits  to  run  from  2  to  12  horizontally  and  from 
3  to  1 1  vertically.  Exit  this  screen. 

7.  Choose  “CONTOUR  LEVELS.”  Change  the  number  of  contours  to  14.  The  contour 
levels  are  recalculated  instantly,  and  the  screen  is  updated  to  show  14  levels. 

Change  contour  levels  13  and  14  to  +0.2,  and  -0.2  volts,  respectively.  Exit  the 
screen. 

8.  Choose  “CONTOUR  MARKS.”  Change  “Number  of  Mark.s”  to  3.  The  screen  is 
redrawn  to  show  3  marked  contours.  Change  the  second  marked  contour  value  to  -.2 
volts  and  the  mark  to  Change  the  third  marked  contour  value  tc  +.2  volts,  and 
the  mark  to  “+”.  Exit  the  screen. 

9.  Choose  “LABELS  AND  INCLUDES.”  Change  the  “Title”  to  “Biasing  Cubes 
Example.”  Exit  the  screen. 

10.  Select  “OTHER  OPTIONS.”  If  you  are  on  a  color  tenninal,  toggle  that  option  to 
“YES.”  Exit  this  screen. 

1 1 .  Return  to  the  pull-down  menu,  and  .select  “Make  Plot”  and  then  “Show  Plot.”  A 
black-and-white  version  of  the  contour  plot  appears  in  figure  5.8.  Hit  <RET''  to  go 
back  to  the  menu.  Exit  “Plot”  menu,  and  then  exit  Scanner. 
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Note  the  prcscnv  >>t  the  tiles  (01^.2’'  and  ■‘2Cubcs.scan.”  The  plot  data  is  contained  in 
“ton. 2'  and  can  fx-  displaced  or  printed  using  any  of  the  DynaPAC  graphical  interfaces. 
“2Cubes.scan  stores  the  options  used  when  running  the  data  scanner,  which  will  be 
restored  in  the  next  Scanner  run  in  this  directory. 

5.9  Running  Particle  (Jenerator  (PartGen) 

Issue  the  command  “DynaPre”  to  generate  input  file  to  module  PartGen.  From  the  main 
menu,  choose  ‘PartGen”  and  then  select  “Edit  Script.”  The  screen  will  appear  as  in  figure 
5.9.  Change  the  “Comment”  to  “Biasing  Cubes  Sample  Problem.”  Toggle  “Particle 
Type”  to  “CONTOUR.”  A  pop-up  form  will  come  up  showing  the  default  cut-plane 
orientation.  Change  “Contour  Potential”  to  “0.1”  volts,  and  change  “Cut  Offset”  to  “7.0” 
(middle  Z-plane).  Exit  this  form.  Now,  change  the  “Number  of  species”  to  “2,”  The 
screen  will  change  instantly  to  show  two  species  (O"^  and  c  as  defaults).  Exit  this  form. 
Select  “Make  Script”  from  the  menu  to  write  out  file  “pg_in.”  Exit  “PartGen”  and  exit 
“DynaPre”.  Inspect  the  file  “pg_in.”  At  the  UNIX  prompt,  type; 

%PartGen_SUN4  <  pg_in  >  pg_out  & 

The  output  file,  “pg_out,"  echoes  the  input  file  “pg_in,”  shows  general  problem  and  grid 
parameters,  and  details  the  particles  generated. 

5.10  Running  Particle  Tracker  (Tracker) 

Run  “DynaPre”  again  to  generate  input  file  to  module  Tracker.  Choose  “Tracker”  from 
the  main  menu  and  then  select  “Edit  Script.”  The  screen  will  appear  a.s  in  figure  5.10. 

Change  the  “Comment”  to  “Biasing  Cubes  Sample  Problem.”  Change  “Title”  to  “Plotting 
Particles  at  sheath=0. 1  volts.”  Change  “Number  of  Species”  to  “2”  to  specify  two 
different  species  (species  1  and  species  2).  Toggle  “Process”  to  “PLOT_PARTICLES.”  A 
new  form  is  popped  up.  Change  “Cut  Offset”  to  “7.0”  ^middle  Z-plane).  Exit  this  form 
The  screen  will  change  instantly  to  show  defaults  corresponding  to  the  selected  option 
"PLOT_PARTICLES.”  Change  “Plot  Limit”  to  run  from  2  to  12  in  X  direction  and  from 
3  to  1 1  in  Y  direction.  The  screen  now  appears  as  in  figure  5.1 1.  Exit  this  form  and  select 
"Make  Script”  to  write  out  file  “tr_plot_in.” 
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Select  “Edit  Script”  again  to  return  to  “Edit"  form.  We  are  going  to  create  yet  another 
input  file  to  Tracker.  Change  “Title”  to  “Tracking  Panicles  from  sheath=0.2  volts.” 
Toggle  “Process”  to  “TRAJECTORIES”  and  accept  the  projection  plane.  The  screen  u  Ul 
change  to  show  defaults  corresponding  to  the  selected  option  "TRAJECTORIES.” 
Change  “B  Field”  to  be  “4.e-5  tesla”  in  “Z”  direction.  Increase  “Max  Steps”  from  500  to 
1000  and  change  “DX_MAX”  from  .3  to  .1  for  better  electron  trajectories.  The  screen 
now  appears  as  in  figure  5.12.  Exit  this  form  and  select  “Make  Script"  to  write  out  file 
“tr^trajjn.”  Exit  “DynaPre.” 

Now,  at  the  UNIX  prompt,  issue  the  command: 

%Tracker_SUN4  <  tr_plot_in  >  tr_plot_out  &  | 

Then  run  your  favorite  “fon.2  plot-reader”  to  display  the  plot  as  in  figure  5. 1 3.  On  color 
plot,  ions  are  shown  in  red  and  electrons  are  shown  in  green.  Note  that  panicles  in  bipolar 
zones  are  not  generated. 

You  might  want  to  save  this  fon.2  file  as  it  will  be  overwritten  by  the  next  command.  At 
the  UNIX  prompt,  type  the  command 

%Tracker_SUN4  <  tr_traj  jn  >  tr__traj_out  & 

Again,  run  your  favorite  “fon.2  plot_reader”  to  display  the  panicle  trajectories  as  in 
figure  5.14. 

Note  that  ion  panicles  got  sucked  right  into  the  object  while  electrons  were  just 
oscillating  about  the  magnetic  field  until  they  fell  into  high-field  region  between  the  two 
cubes  where  they  were  pulled  into  the  cube’s  face. 

Some  new  subfiles  have  been  created  by  PanGen  and  Tracker  here: 
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2Cubes.Fn 

2Cubes.P0TG 

2Cubes.POTS 

2Cubes.CUR 


Contains  particle  atinbutes,  created  by  PanGen  and 
updated  by  Tracker  as  panicles  are  tracked. 

Contains  time-dependent  spatial  potentials. 

Contains  time-dependent  surface  potentials. 

Contains  time-dependent  current  to  surface  information. 


It  is  interesting  to  look  at  trajectories  of  fewer  electron  particles  instead.  Now,  rerun 
“PanGen”  with  command: 


%PanGen_SCN4  <  pg_in  >  pg,_out 

Although  you  could  run  “DynaPre”  again  to  build  a  new  “tr_traj_in’',  it  is  simple  enough 
to  make  the  change  manually  using  a  text  editor.  Edit  the  following  lines  in  “tr_traj_in”  to 
read: 


XJJMIT  7.0  7.5 
Y_LIMIT  7.0  1 3.0 
Z_LIMIT  0.0  0.0 

Then,  run  “Tracker”  again  with  command: 

‘7cTracker_SUN4  <  tr  _traj..in  >  tr_traj_out 

Once  more,  run  your  favorite  plotter  to  display  the  new’  trajectories,  which  will  be  as  in 
figure  5.15 

Now,  just  for  fun,  let’s  change  the  magnetic  field  to  be  along  X-axis  and  track  sheath 
panicles  on  the  X=8  plane  and  within  5.5  <  Y  <  6.  5  <  Z  <  6  range.  Run  “DynaPre”  to 
build  appropriate  “pg_in”  and  “tr  _trdj_in”  input  files  to  “PanGen”  and  "Tracker,” 
respectively.  Then  run  “PanGen”  and  “Tracker”,  e.j^. 

%PartGen_SUN4  <  pg..in  >  pg_out 
%Tracker  ,SUN4  <  tr_traj_in  >  tr_traj_out 

finally,  plot  the  trajectories,  which  will  look  like  figure  5.16. 
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5.11  Summary 


This  example  has  taken  you  through  most  of  DynaPAC  modules;  defining  an  object; 
setting  up  computational  grids;  calculating  and  displaying  potentials  in  the  surrounding 
space;  generating,  tracking  particles,  and  displaying  particle’s  trajectories. 


dynapac 


Exit  Fil=a 


utter  Show  Write  Beset 


.-.e^p 


Modifying  PRIMABV  or  id 
Or  id  dimension:  m  s  5 

Mash  size:  Q.200Q 

Object  extents  in  grid  units: 

1.2500  <-  X  <-  4.7500 

2.2500  <-  Y  <-  3.7500 

2.2500  <-  2  <-  3.7500 


tR.ii^NSLATE  OBJECT 
PLOT 


ACCEPT  CHAJIGE3 
IGNORE  CHANGES 
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Figure  5.3  GridTool  plot  for  the  2-Cubes  sample  problem 
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Figure  5.7  Scanner  screen  for  the  Plot/Edit  Plot  option. 
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Figure  5.9  DynaPre/PanGen  editing  screen. 
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Figure  5.14  Particle  trajectories. 
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Figure  5.15  Electron  trajectones. 


108 


6.  Docunientatiun  for  Two-D  Klectroslalic  Code  "(Jilberf 


6.1.  Introduction 

Gilheri  is  a  general,  two-dimensional  (R-Z  or  X-Yi  tinile-elcment  electrostatic 
code.  It  can  solve  a  fairly  general  version  of  Foiss^ai’s  equation,  and  can  emit  and 
track  charged  particles  in  electric  and  magnetic  tields. 

Poi,sson’s  equation  is  solved  taking  into  account  varving  dielectric  constants  in 
different  regions.  Roundars-  conditions  include  li.xed  potential  conductors,  floating  con¬ 
ductors,  and  intcmal  "w'indowscreen''  boundary  conditions.  .Space  charge  is  contri¬ 
buted  by  charged  particles  iKing  tracked,  by  charged  panicles  deposited  on  insulating 
suifaces,  by  conductivity  effects,  and  by  analytic  "plasma"  formulas. 

Charged  panicles  may  be  normally  emitted  from  conducting  surfaces  ai  fixed 
energy  and  current  density.  Zero  energy  panicles  will  be  emitted  in  space-charge- 
limited  fashion.  An  applied  magnetic  lield  in  the  Z  direction  (tor  eitiicr  geometry') 
may  be  specified,  lip  to  five  dillereiu  particle  species  may  be  specified.  Initial  pani¬ 
cles  may  be.  generated  using  the  auxiliary  routine  ItiiiP 

All  infomiation  concerning  gnddmg.  particle  properties,  input  variables,  potential 
values,  etc.  is  maintained  in  a  single  database  (stored  as  !bn.2()).  Auxiliary  programs 
to  make  plots  or  to  otherwise  query  the  database  are  included  in  the  Gilbert  package. 
Additional  special-purjiose  programs  may  easily  be  written  by  users. 

Presently  Gilbert  is  implemented  under  UNIX  on  a  Ridge  320()  (space).  How¬ 
ever,  it  should  be  readily  convenible  to  any  other  UNIX  machine,  or  to  VAX/VM,S. 
Graphical  programs  are  at  present  designed  for  PostvSeript  printers,  although  they  can 
also  be  di.splayed  on  Tektronix  4014  terminals  or  emulators.  Color  output  is  aniici- 

patpd. 

6.2.  Grid  Generation 

6.2.1.  (iHdding  Requirements 

Gilbert  operates  on  a  finite  element  grid  composed  ol  an  arbitrary  combination  of 
linear  triangles,  bilinear  (4-node)  quadrilaterals,  and  biquadratic  (8-node)  quadrilaterals. 
(In  general,  the  8-n<x!e  quads  will  give  superior  results.)  These  polygons  are  referred  to 
as  "elements”,  and  their  comers  (or  mid-side  points  for  the  8-node  quads)  are  referred 
to  as  "nodes”,  liach  element  lias  assiK’iated  with  it  an  integer  propeily  (KNkl.S) 
called  a  material  number.  Material  1  is  considered  free  space,  and  supports  prop.iga 
tion  of  charged  particles.  Materials  2  1.3  are  considered  dielectrics,  and  are  capable  of 
charge  deposition  on  their  surfaces.  Uach  node  has  as.socKited  with  it  an  integer  pro¬ 
perty  (()<B<2())  called  a  "txiundai^  condition".  Ntxlcs  with  kiundary  condition  0  are 
frce-space  nodes.  Nixies  with  Ixnindary'  conditions  1-20  are  conductor  nodes,  with 
user-specified  jiropertics  (see  below).  [Tcc-space  nodes  on  the  boundary  of  the  grid 


will  have  zero-normal-electnc-ficld  boundary  condiuons. 

An  external  mesh  generator  (such  as  Fatran)  is  used  to  generate  the  grid.  The 
mesh  generator  must  provide  a  sequential  list  of  nodes,  with  their  Iwations  and  boun¬ 
dary'  conditions,  and  a  sequential  list  of  elements,  with  their  detining  nodes  and 
material  numbers.  Elements  must  be  compatible,  i.e..  a  side  of  a  polygon  (or  half-side 
of  an  8-node  quad)  must  also  be  a  side  of  another  polygon  or  it  u  ill  be  considered  a 
boundary.  (Note  that  in  a  triangular  region  some  mesh  generators  will  supply  6-node 
triangles  where  8-node  quads  have  been  requested.  Such  elements  are  not  currently 
acceptable  to  Gilbert,  and  must  be  rezoned  by  the  user.)  The  mesh  generator  should  be 
used  to  eliminate  unintended  internal  boundaries.  AKso,  the  mesh  generator  should  be 
used  to  minimize  abrupt  changes  in  grid  spacing,  and  to  eliminate  particularly  badly 
misshapen  elements. 

6.2.2.  Grid  Generation  using  Patran 

Pending  further  request,  the  only  available  mesh  generator  interface,  ReadPai. 
reads  the  Patran  "neutral  lile".  Patran  (a  prcxluct  of  PDA  Engineering)  is  available  at 
S-CL'BED  on  the  Silicon  Graphics  Iris,  and  may  be  run  with  graphical  interaction  on 
the  Iris  console,  or  in  "BAT"  mode  from  any  other  terminal.  We  also  have  .software 
capable  of  converting  an  IGIiS  file  to  a  Patran  neutral  file.  Consult  Myron  Mandell  or 
Victoria  Davis  concerning  access  to  Patran  or  Patran  documentation. 

For  Gilbert  we  use  Patran  to  construct  a  two-dimensional  mesh  in  the  X-Y  plane 
(i.e.,  the  default  plane  of  the  screen).  The  nodes  may  later  be  .scaled  (to  account  for 
units),  and  interpreted  as  X-Y,  Y-X,  R-Z,  or  Z-R.  Element  material  numbers 
correspond  the  Patran  .MIDs.  and  node  conductor  numbers  correspond  to  the  Patran 
node  configuration  number. 

It  is  not  our  intent  here  to  provide  a  Patran  tutorial,  but  we  will  give  a  brief 
description  of  a  few  of  the  concepts  and  commands  which  might  be  used  to  generate  a 
grid. 

(Brief  Patran  Discussion  to  be  inserted) 

It  is  not  necessary  to  assure  that  all  elements  are  counterclockwise,  as  this  opera¬ 
tion  will  be  performed  by  ReadPat. 

6.2.3.  Grid  PrtK:es.sing  wilh  ReadPat 

ReadPat  creates  and  initializes  the  database  (fort. 20)  with  the  gridding  infonna- 
tion  and  with  default  values  for  vtirious  other  options  and  parameters.  ReadPat 
expects  to  find  the  Patran  neutral  hie  on  logical  unit  8  (fon.8).  It  also  accepts  certain 
options  from  standard  input,  and  wntes  an  output  file  on  standard  output. 

The  options  accepted  by  ReadPat  on  standard  input  are; 
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Input  Syntax 

Meaning 

Default 

DPRN 

Diagnostic  Print 

Omit  Diagnostic  lYint 

ECHO 

Echo  Patran  Neutral  File 

No  Echo 

NOEC 

Do  not  echo  Patran 

Neutral  File 

No  Echo 

FLIP 

Interpret  input  as 

Interpret  input  as 

Y-X  or  7.  R 

X-Y  or  R-Z 

UNIT  r 

Input  unit  is  r  meters 

r  -  1 . 

UNIT  INCH 

Input  unit  is  .0254  meters 

UNIT  FEE'r 

Input  unit  is  ..^048  meters 

UNi  r  FOOT 

Input  unit  is  .3048  meters 

UNIT  CENT 

Input  unit  is  .0100  meters 

UNIT  METE 

Input  unit  is  1  meter 

UNIT  MILL 

Input  unit  is  .(X)l0  meters 

DIFL  i  r 

Material  i  has 

dielectric  constant  r 

Dielectric  constant  -  ' 

SIGM  i  r 

Material  i  has 
conductivity  r 

Conductivity  =  0. 

RZ 

(No  Effect) 

Geometry  is  R-Z 

XY 

Geometry  is  X-Y 

unless  Rmin  <  0 

END 

No  further  input 

6.2.4.  Grid  Illustration  with  GridrUn 

After  creating  the  database  with  ReadPat,  it  is  possible  to  produce  PostScript  plots 
of  the  interpreted  grid  using  the  GridPloi  program.  GridPloi  will  read  the  database 
(fort.20)  and  create  graphics  commands  on  logical  unit  2  (fort. 2).  PonPloi  will 
translate  these  commands  into  PostScript  commands  on  its  standard  output,  which  may 
be  piped  or  otherwise  sent  to  a  convenient  PostScript  printer.  Alternatively,  Tek40N 
will  write  Tektroni.\  4014  commands  on  its  standard  output. 

GridPloi  will  draw  element  outlines  and  node  points,  and,  for  R-Z  geometry,  will 
draw  the  axis  of  symmetry.  Elements  of  materials  2-15  will  be  shaded  in  gray.  At 
user  preference,  GridPloi  can  print  element  numbers,  element  material  numbers,  con¬ 
ductor  numbers,  or  node  numbers.  To  accept  a  choice,  respond  to  the  prompt  with  a 
simple  carriage  return  <CR>.  Any  non-null  response  will  reject  the  choice. 

GridPloi  has  a  windowing  capability,  allowing  plots  to  be  restricted  to  only  a  por¬ 
tion  of  the  grid. 
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6.3.  Running  Gilbert 


6.3.1.  General  Considerations 

Gilbert  is  run  in  the  directory  where  its  database  (fon,20)  resides.  It  reads  an 
input  tile  from  standard  input,  writes  output  on  standard  output,  ami  continually 
updates  its  database.  Calculation  takes  place  in  four  phases;  { 1 )  read  input  and  related 
initialization;  (2)  generate  Laplace  equation  matrices  (if  necessary);  (3)  initialize  poten¬ 
tials  (if  necessary);  (4)  perform  iterative  cycles  to  resolve  non-hnear  and/or  time- 
dependent  phenomena. 

Gilbert  first  reads  in  all  current  values  of  options  and  parameters  from  the  data¬ 
base  (fort.20).  It  is  envisioned  that  a  screen-oriented  input  program  will  be  made 
available  to  manipulate  these  options  and  parameters.  At  present,  however,  user- 
modifiable  options  and  parameters  may  be  changed  via  command  from  siandiu-d  input. 
These  options  and  piu'ameters  and  their  input  syntax  are  described  in  the  next  section. 

After  insening  the  mesh  into  the  database  with  ReadPut,  it  is  necessary  to  gen¬ 
erate  the  matrices  to  be  used  to  solve  for  the  potential,  as  well  as  to  perform  certain 
other  initializations.  This  is  done  in  response  to  the  "GenMat"  command.  ("GenMat" 
is  set  to  .True,  by  RcadPat.  It  is  automatically  set  to  .False,  after  the  initialization  is 
perfonned.) 

Potentials  in  space,  and  charges  on  floating  conductors,  are  initialized  in  response 
to  die  command  "IniPot".  The  variable  ''fniPoi"  may  be  set  to  .True,  from  the  input 
stream.  It  will  be  set  to  .False,  following  the  initialization. 

The  variable  "Cycle"  determines  whether  iterative  cycles  are  to  be  perfonned.  it 
is  set  to  .True,  when  a  non-zero  number  of  cycles  is  requested  via  the  input  variable 
"NCYC".  Each  cycle  consists  of  four  phases:  (1)  emit  particles,  track  the  particles, 
and  deposit  charge  accordingly  on  conductors  and  insulating  surfaces;  (2)  computed 
charge  accumulation  due  to  conductivity  effects;  (3)  calculate  space  charge  due  to  par¬ 
ticles,  particle  deposition,  and  plasma  effects,  and  update  the  potential  solution  accord¬ 
ingly;  (4)  generate  written  and  printed  output  as  requested.  If  further  cycles  are 
required,  subsequent  runs  may  resume  from  the  current  state  of  the  database. 


6.3.2.  Common  Block/lnpul  Specifications 
6.3.2. 1.  Boundary  Conditions 

COMMON/BCOND/BCTyp(NCMX),VOLTS(NCMX),  BCPann(NCMX).  Capaci(3.NCMX) 

NCMX  is  a  parameter  denoting  the  maximum  number  of  conductors  (currently 
20).  BClTyp  describes  the  type  of  boundary  condition.  Currently  implemented  values 
are  "FIX  "  for  fixed  potential  conductors,  "FLOQ"  for  charge- initialized  floating 
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conductors,  "FLOV"  for  potentiui-iniiialized  floating  conductors,  "lONO"  to  indicate 
tliat  this  collection  of  nodes  is  not  to  be  treated  as  a  conductor,  and  "SCRFi"  for  win- 
dowscrcen  internal  boundary.  VOl.'fS  describes  the  voltage  of  a  fixed-potential  boun¬ 
dary,  the  initial  voltage  of  a  potential-initialized  floating  Ixiundary,  or  the  nominal  vol¬ 
tage  of  a  windowscreen  boundary.  BCParni  may  contain  either  the  parameter  relating 
effective  voltage  to  field  discontinuity  for  a  windowscreen  boundary,  or  the  initial 
charge  for  a  charge-initialized  floating  boundary.  Capaci  is  an  array  containing  intcr- 
conductor  capacitances.  The  first  two  entries  of  each  triplet  are  integers  giving  the 
conductor  numbers  to  be  capaciiiveiy  connected,  and  the  third  is  the  capacitance  in 
farads  (R-Z  Geometry)  or  farads/meier  (X-Y  Geometry).  Note  that  changes  in  capaci¬ 
tance  values  will  be  implemented  only  during  the  matrix  generation  (GenMatj  phase  of 
the  code. 


Relevant  input  variables 

are: 

Input  Syntax 

Meaning 

De  fault 

FIX  i 

Fixed  potential  on 
conductor  i 

Floating  Potential 

VOLT  i  volts 

Define  potential 
on  conductor  i 

Volts(I)  =  0. 

FLOA  i 

Potential-initialized 
floating  potential 
on  conductor  i 

(Default) 

IGNO  i 

Ignore  specification 
for  conductor  i 

CFIAR  i  coul 

Charge  initialized  floating 
conductor  i,  with  coul 
coulombs  of  charge 

Voltage-initialized 

SCRE  i  alpha 

Windowscreen  internal 

boundary'  on  conductor  i, 
with  parameter  alpha 

Voltage-initialized 

CAPA  id  ic2  ftu'ads 

Interconductor  Capacitance 

Stray  Capac.  only 

CAPA  RESET 

zero  all  capacitances 

6.3.2.2,  Plasma  Parameters 

Common  /Plasma/  Rho,  TO,  RLamda,  Dsq,  Linear.  Planar,  Convrg,F.lec,Ions 
Logical  Linetir,  Planar,  Convrg,Hiec,lons 

This  common  block  contains  the  plasma  density,  temperature,  and  Debye  length. 
Dsq  is  the  square  of  the  Debye  length.  The  three  logical  variables  describe  the  formu¬ 
lation  to  be  used  for  pla.sma  space  charge,  with  the  options  "Linear"  for  linear  (Debye) 
screening,  "Planar"  for  Child-Langmuir  .screening,  "Convrg"  for  spherically  convergent 


(Langmuir-Blodgett)  screening,  "Elec"  for  screening  by  only  the  electron  component  of 

the  plasma  (e.g.,  if  ions 

are  tracked),  and  "Ions" 

for  screening  by  only  the  ion  corn- 

ponent  of  the  plasma  (c.g 

if  electrons  are  tracked) 

.  Relevant  input  variables  are: 

Input  Syntax 

Meaning 

Default 

PLAS  LINE 

Linear  Screening 

No  Screening  or  As  Before 

PLAS  PLAN 

Planar  Screening 

No  Screening  or  As  Before 

PLAS  CONV 

Convergent  Screening 

No  Screening  or  As  Before 

PLAS  ELEC 

Electron  Screening 

No  Screening  or  As  Before 

PLAS  IONS 

Ion  Screening 

No  Screening  or  As  Before 

PLAS  other 

Turn  off  Screening 

As  Before 

DEBY  length 

Set  Debye  Length; 

7430  meters,  or  As  Before 

Redefine  RHO 

1  m  ■ ,  or  As  Before 

RHO  density 

Set  RHO; 

Redefine  RLamda 

TEMP  temp 

Set  TO; 

1  eV,  or  As  Before 

Redefine  RLamda 

6. 3.2.3.  Mesh  Parameters 

COM.MON/.MESH/NNc)dcs,NHIt,IB.Set,R.Mm.RMax,ZMin,Z.Max.RZGcom,XYGeom 
Logical  RZGeom.XYGcom 

NNodes  and  NElt  describe  the  number  of  nodes  (maximum  3100)  and  elements 
(maximum  3000)  actually  in  the  grid.  RMin,  RMax,  ZMin,  and  Z.Max  describe  the 
minimum  and  maximum  R  (or  X)  and  Z  (or  Y)  coordinates  of  the  nodes.  The  logical 
\ariables  RZGeom  and  XYGeom  describe  whether  the  geometry  is  to  be  treated  as  R- 
Z  or  X-Y.  (Note  that  this  common  block  contains  substantial  further  gridding  infoimu- 
non,  which  is  omitted  here.) 

Relevant  input  keywords  are  "XY"  and  "RZ",  which  reset  the  geomctiw  type. 
(Note  that  if  the  geometry  type  is  actually  reset,  the  matrix  initialization  ("GenMat") 
sihould  be  redone.) 

6. 3. 2. 4.  Particle  Parameters,  Magnetic  Field,  Time.step 

CO.MMON  /P.ARCON/  EQ(f>h  H.M(5).  F.By.Vl(.3j,  BZ,  Omega(S),  DeiT,  MiiTor 
Logical  Mirror 

EQ  and  EM  give  the  (signed)  particle  charge  (coulombs)  and  mass  (kilograms) 
L)r  tlie  .species  of  particles  of  each  species  to  be  tracked.  EByM  -  f3Q/EM.  BZ  is  the 
z  componeni  of  magnetic  field  (W-m  ").  Omega  -  EByM*BZ.  DelT  is  the  timestep 
()cr  cycle,  and  also  the  ma.ximum  timestep  per  particle  push.  Mirror  indicates  that  Z=() 
(for  R  Z)  or  Y=0  (for  X-Y)  is  a  mirror  plane.  Relevant  input  variables  (applied  to  the 
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current  particle  species  number)  are: 


Input  Syntax 
PART  SPEC  n 

PART  CHAR  coul 
PART  MASS  kg 
PART  MASS  a  AMU 
BZ  w 

BZ  g  GAUS 
DELT  sec 
MIRROR 
MIRROR  OFF 


Meaning 
Particle  Species 

lor  subsequent  input 
Panicle  charge  (couioinbs) 
Particle  mass  (kilograms) 
Panicle  mass  (amu) 

.9 

Magnetic  field  (W-m  ") 
Magnetic  field  (gauss) 
Timestep  (seconds) 
Declare  Mirror  Plane 
Omit  Mirror  Plane 


Default 

1 

-i.bxio;^*^ 
9.1x10'"’*  kg 
9.1x10'-^*  kg 
0. 

0. 

-O 

l.xlO  ^ 

No  Mirror 
No  Mirror 


6.3.2.S.  Emission  parameters 

Common  /EMISSN/  Window(4),CurDen(NCMX),Energy(NCMX),NPSEG(NCMX), 

*  JSpeces(NCMX),NPMod(NCMX) 

Window  give  the  XMin,  XMax,  YMin,  YMax  from  which  emission  can  take 
place.  (This  is  useful  if  only  a  portion  of  a  conductor  is  emitting.)  CurDen  gives  the 
maximum  current  density  (A-m  *■)  emitted  from  each  conductor.  Energy  gives  the 
energy  at  which  panicles  are  emitted,  with  zero  energy  indicating  space-charge-limited 
emission.  NPSeg  gives  the  number  of  panicles  to  .9e  emiiied  for  each  line -segment  of 
conductor.  JSpeces  gives  the  .species  nun;ber  of  particles  to  be  ontitted  from  each  con¬ 
ductor,  as  specified  by  the  most  recent  "PART  SPEC  n"  card.  (See  panicle 
specifications  section  above.)  NPMod  gives  the  number  of  timesteps  between  emis¬ 
sions.  (Useful  for  problems  having  both  slow  and  fast  panicles.) 

Relevant  input  ptuameters  are: 

Input  Syntax  Meaning  Dclauli 

EMIS  WIND  xl  x2  yl  y2  Delinc  emission  window  Entire  Gnd 

EMIT  COND  i  ENER  c  CURD  j  NPSE  n  Mod  m  Define  emission  (No  emission) 

propcnics  for 
conductor  i 


6.3.2.6.  Caleuiation  Directives 

Common  AVhToDo/  GenMat,  IniPot,  Cycle, LstCyc,  NCyc,  Time 
Logical  GenMat,  IniI\M,  Cycle 

GenMat  defines  whether  Laplace  matrix  initialization  is  to  be  performed.  IniPot 
defines  whether  potential  initialization  is  to  be  performed.  Cycle  defines  whether  itera¬ 
tive  cycles  are  to  be  performed.  LstCyc  defines  the  number  of  the  previous  cycle  or 
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tiniesiep.  NCyc  definc.s  ilic  number  of  cycles  to  tx;  performed,  'l  ime  is  the  running 
total  of  the  DelT’s  for  each  iimestep.  Input  parameters  :ue; 


input  Syntax 

Meaning 

Dclault 

NCYC  n 

Number  of  cycles; 

Set  CYCLE  to  .True. 

.No  cycles 

GENM 

Initialize  Laplace  matrix 

Set  to  .  rrue  by  ReadPai 
Set  to  -f  alse,  after  done 

INiP 

Initialize  Potentials 

.False. 

CYCLE  RESET 

Set  LstCye='rime=() 
Maintain  existing  panicles 
and  charge  deposition. 

Continue  from  previous 

6. 3.2.7.  Particle  Tracking  Parameters 

Common  mack/  XYNoB,  XYB.RZ.NoB.  RZB.DXMax,DYMax,DF;Ma.v 
Logical  XYNoB, XYB.RZNoB.RZB 


'I  he  four  logical  variables  are  set  iniernaily  to  determine  the  appropriate  type  of 
panicle  pushing.  DXMa.x  and  DYMax  define  I’.v  mavimum  X  (or  R)  and  Y  (or  Z) 
distance  a  panicle  can  move  in  a  om.  ,i^p.  If  this  is  exceeded  with  timestep  DelT, 
then  multiple  particle  timesteps  will  be  taken  within  the  cycle  for  that  panicle.  Simi¬ 
larly,  DEMax  is  the  maxitmin  kinetic  energy  change  per  particle  push.  Relevant  input 
v.iriabies  are; 


Input  Syntax 
DXMA  dx 
DY.MA  dy 
DEMA  de 


Meaning 

Maximum  DX  per  panicle  push 
.Maximum  DY  per  panicle  push 
Maximum  kinetic  energy 
change  per  panicle  push 


Default 

O.lv(RMax-RMin) 
0. 1  <(ZMax-ZMin) 
1x10^  eV 


6. 3.2.8.  Other  Miscellaneous  parameters 

Common  /MISC/  Maxitr,  IPrPar,  ParPsh,  Conduc 
Logical  ParPsh,  Conduc 

Maxitr  is  the  maximum  number  of  iterations  within  the  ICCG  potential  solver. 
IPrPar  turns  on  panicle  diagnostics.  ParPsh  defines  whether  any  panicle  pushing  is 
required.  (ParPsh  is  set  to  .'I’rue.  when  emission  is  requested,  but  may  also  be  set  by 
the  user.)  Conduc  defines  whether  conductivity  is  to  be  considered.  (Conduc  is  ini¬ 
tially  set  to  .False,  by  RcadPai,  and  changed  to  .True,  wl.en  ReadPat  encounters  a 
’SIGM’  input  card  assigning  conductivity  to  a  material.  Relevant  input  variables  are: 

Input  Syntax  Meaning  Default 
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MAXI  II 
IPKP  II 
PARF 
PARP  OFF 
COND 
COND  OFF 


Maxiiiuim  number  ot  lC'C(i  iicraiioii^ 

Piirlii  le  Diagnoslics 

Turn  on  particle  pushing 

I'lirn  olT  panicle  pushing 

Turn  on  conductivity 

Turn  off  conductivity 


30 

0 


6.3.2, 9.  O’ttput  Specifications 

It  is  commonly  desired  to  print  or  plot  a  quantity  as  a  function  of  time  or  cycle 
number.  When  Subroutine  Input  encounters  an  ’OUTP'  card,  it  transfers  control  to 
Subroutine  OutSpc  to  interpret  the  remainder  of  the  card.  The  capabilities  and  syn¬ 
taxes  of  this  feature  are  discussed  below  in  the  section  titled  Output  Discussion  . 


6.3.3.  Windowscreen  Boundary  Condition 

We  have  incorporated  into  tlic  code  an  algorithm  for  treating  the  w indowscreen  as 
an  interior  boundary  condition.  The  inciluKl  dctennincs  the  IikuI  elfcctivc  potential  ol 
the  windowscreen  sclf-eonsistemly  with  the  calculated  electric  fields  on  either  side  of 
the  .screen. 

Normally,  the  windowscreen  spacing  is  smaller  than  any  other  parameter  (resolu¬ 
tion  or  Debye  Length)  in  the  problem.  In  this  case,  the  difference  between  the 
effective  potential  of  the  windowscreen  and  its  actual  (circuit)  potential  is  proponional 
to  the  discontinuity  of  electric  field  across  the  screen; 


V^f^.-V^  =  -a(T)L(E^^-E,).n 

where  L  is  the  interwire  spacing,  af  T)  is  a  function  of  the  transparency  of  the  screen, 
b  above  and  l)elow  the  screen,  and  n  is  the  upward  unit  nor¬ 

mal  to  the  screen.  The  ptu'ameter  ("alpha")  which  must  be  entered  when  specifying 


the  ’’SCRE"  boundary'  condition  is  a('l')AL. 
various  transparencies. 


The  table  below  gives  values  of  a(T)  for 


Transjiarency 

a(T) 

5m 

.01X3 

60% 

.0310 

70% 

.04X8 

«()% 

.0734 

m 

.120 

95% 

.174 
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6.3.4.  L'ser-Modifiable  Kuutine.s 


It  is  not  possible  to  oiilcuiatc  in  n  general  way  every  variation  on  ihc  eleelrosiaiic 
theme.  However,  we  have  med  to  amicipaie  some  of  the  areas  in  which  users  wall 
wish  to  "hardwire"  ptirameter.s  and  functional  fonns  for  their  particular  applications 
and  to  localize  these  in  user  modifiable  routines.  'I'he  routines  intended  for  such  pur 
poses  are  extensively  internally  documented.  We  list  these  routines  below; 

6. 3.4.1.  Subroutine  ExtCur 

ExtCur  provides  the  capability  of  imposing  a  time-dependent  potential  on  "MX  " 
and  "SCRE"  conductors,  or  providing  an  external  current  (which  may  depend  on  time 
and  potential)  to  floating  conductors. 

6.3.4.2.  Subroutine  Output 

Output  provides  the  capability  of  printing  or  writing  for  later  processing  several 
categories  of  quantities  in  a  general,  automated  way.  However,  the  user  may  have 
specific,  unanticipated  quantities  of  which  he  would  like  to  keep  track.  A  stencil  is 
provided  in  Subroutine  Output  for  the  user  to  calculate  and  pnnt/writc  such  a  quantity. 

6,3.5.  Output  Discussion 

6.3.5. 1.  User  Requested  Output 

The  user  may  request  output  to  be  printed  or  to  be  w'ritten  on  itnfomiaited  Hie 
fort. 30.  The  output  nia)  bo  conductor  or  node  charges,  conductor  or  node  potentials, 
nodal  fields,  and  total  panicle  charge  or  dipole  moment.  A  program  LinFl:  is  provided 
to  plot  the  data  on  fort. 30,  As  with  other  computational  parameters,  the  output 
requests  are  stored  in  common  blocks  and  maintained  in  the  database  for  restan  pur¬ 
poses.  We  will,  as  done  with  input  parameters  above,  categorize  the  output  requests 
by  common  block. 

Changes  in  output  requests  should  be  done  only  when  the  INIP  request  is  used  to 
establish  the  calculation  at  cycle  1. 


6.3.5. 1.1.  (icneral  Requests 

Common  /OMods/  ModPn,  ModWn 


The  cards 


OUTPUT  MODP  npriiit 
Oin  PUT  MODW  nwnte 

specify  that  pnnting/wriiing  takes  place  every  nprim/nwriie  cycles.  A  value  of  zero 
indicates  no  printing/writing,  'llie  default  is  to  set  btith  values  to  zero.  Al.so,  the  card 


UK 


OU  i’PUT  RESii'r 


sets  the  default  values  for  all  output  parameters  except  ModFn  and  ModWn, 

6.3.5. 1.2.  Conductor-Related  Quantities 

Common  /OCnd/  lConds(2{)),  l.Chrge,  Lf^ot,  Lr.mit,  LRetrn 
Logical  l.Chrge,  I. Pot,  IJimii,  LRetrn 

The  values  in  the  above  common  block  are  controlled  by  the  cards 

OUTPUT  CONDUCrOR  icond  lOFI'i 
OUTPUT  CONDl'CrOR  CHARGE  !OFL| 

OUTPUT  CONDUCTOR  POTENTIAL.  lOEl'l 
OUTPUT  CONDUCTOR  EMISSION  (OFf  J 
OUTPUT  CONDUCTOR  RETURN  |OFF| 

The  default  is  to  phnt/write  potential,  not  to  pnnt/wriic  charge,  cumulative  emission 
current,  or  cumulative  return  current,  and  to  have  no  conductors  specified. 

6.3.5. 1.3.  Nodal  Potentials 

Common  /OPots/  .NodPoi{20) 

Nodal  potentials  are  requested  by  the  card 

OUTPUT  POTFIN  TIAI.  NODE  inode  |Oi  11 
By  default,  no  nodes  are  specilied. 

6.3.5. 1.4.  Nodal  Fields 

Common  /OFlds/  NodFld(2,20),  LEI.  LE2,  LETot 
Logical  LEI,  LE2.  LETot 

Nodal  fields  are  requested  by  the  cards 

OUTPUT  FIELD  NODE  inode  fELEMENT  ieltl 
OUTPUT  FIELD  NODE  inode  [OFF! 

OUTPUT  FIELD  lOTAL  [OFF] 

OUTPUT  FIELD  dir  [OFF] 

where  dir  is  one  of  |R,  X,  1  ]  to  prini/write  the  R  or  X  field  component,  and  one  of  [Z, 
Y,  2]  to  print/write  the  Z  or  Y  field  component.  If  no  element  i.-*  specified  for  a  node, 
Gilbert  will  deiennine  an  appropriate  element.  The  default  is  to  have  no  nodes  or  de¬ 
ments  specified,  and  to  print/write  only  the  total  electric  field. 
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6.3.5.I.5.  Particle  Quantities 


Common  /OPurls/  LQPan(0.t,  Li)ip2(0) 

Logical  LQPan,  LDipl.  LDip2 

Particle  quantities  itfe  requested  by  the  cards 

OL  TPUT  PARTICLL  SPHC'ILS  n 
OUTPUT  PARTICLL  C'UARGL  10LL| 

OU  I  PUT  PARTICl.K  DIPOl.L  dir  lOI  Ll 

where  a  value  n  from  1  to  5  specilies  a  particle  species  for  subsequent  OUTPUT  PAR¬ 
TICLE  commands  (any  other  value  indicating  sum  over  all  particles),  and  dir  is  one  of 
[R,  X,  Ij  to  print/write  the  R  or  X  dipole  component,  and  one  of  [Z,  Y,  2|  to 
print/write  the  Z  or  Y  dipole  component.  The  default  is  to  omit  calculation  ol  particle 
quantities. 

6.4.  Auxiliary  Programs 

6.4.1.  GridPlot 

GridPlot  has  been  described  above. 

6.4.2.  InitP 

InilP  allows  the  user  to  specify  initial  particles  throughout  the  grid.  I'he  u.ser  will 
be  prompted  for  a  density  for  each  particle  species.  If  a  positive  density  is  given,  Iniil^ 
will  place  one  macroparticle  on  each  triangle  or  4-node  quadrilateral,  and  four  macro¬ 
particles  on  each  8-node  quad.  InitP  should  logically  be  mn  following  an  initial 
potential  (I.NIP)  calculation  by  Gilbert. 

6.4.3.  Contours 

ConUturs  allows  the  user  to  plot  equipotential  contours  in  the  grid.  Contour 
v.ilucs  may  be  either  quasi-logariihmically  spaced  (1,  2,  .5,  10,  20,  ...)  or  evenly 
spaced.  Windowing  capability  allows  plotting  to  be  limited  to  a  portion  of  the  grid. 

6.4.4.  Scatter 

Scatter  plot.s  the  positions  of  particles  ul  each  species  currently  in  the  problem. 

6.4.5.  LinPIt 

LinPli  allows  the  user  to  plot  quantities  written  to  tile  fon.,M)  via  the  output 
options.  The  quantities  may  be  plotted  either  vs.  time  or  vs.  cycle  number,  and  a 
range  of  cycles  may  be  selected.  Also,  the  user  may  select  the  quantities  to  be  plotted. 
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The  user  may  readily  create  custom  versions  of  LuiFli  to  plot  multiple  curves  on 
a  plot  or  to  plot  derived  quantities. 


1 
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